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1.  Introduction 


This  report  describes  the  development  of  a  semi-active  laser  (SAL)  terminal  guidance  model. 
SAL  guidance  typically  consists  of  a  scout  illuminating  a  target  with  short,  high-energy  laser 
pulses  in  a  near-infrared  (IR)  wavelength  (figure  1).  Some  of  the  energy  from  each  pulse  reflects 
off  of  the  target  and  into  the  seeker’s  multi-section  IR  detector.  By  comparing  the  power 
measured  in  each  section  of  the  IR  detector,  the  seeker  approximates  the  location  of  the  reflected 
power  within  its  field  of  view  (FOV).*  Using  this  information,  the  seeker  returns  control  signals 
to  the  projectile’s  maneuver  system  to  steer  towards  the  target. 


In  this  model,  the  seeker  is  assumed  to  be  inside  the  nose  of  the  projectile. 
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The  first  objective  of  the  model  is  to  determine  the  laser  power  distribution  at  the  seeker.  This  is 
accomplished  by  calculating  the  laser  beam’s  transmission  path  and  power  loss.  The  beam  is 
modeled  stochastically  through  a  set  of  rays  that  form  a  solid  cone  of  the  given  divergence.  Each 
ray’s  transmission  path  and  power  loss  is  calculated,  and  its  power  is  summed  with  the  other 
rays’  contributions  at  the  seeker.  Ray  tracing  determines  ray  transmission  paths  from  designator 
to  target,  and  target  to  seeker.  The  Beer-Lambert  Law  is  used  to  compute  power  loss  due  to 
atmospheric  attenuation.  Lambert’s  Cosine  Law  and  designator-seeker-target  geometry  are  used 
to  determine  the  portion  of  the  power  that  the  seeker  receives  after  it  reflects  off  of  the  target. 

The  second  objective  of  the  model  is  to  convert  the  seeker  power  distribution  into  projectile 
guidance  signals.  The  guidance  signals  are  calculated  by  dividing  the  power  received  in  the 
pitching  and  yawing  directions  by  the  total  power  received  by  the  seeker.  The  guidance  signals 
are  not  considered  reliable  unless  the  seeker  power  is  above  a  signal-to-noise  (S/N)  threshold. 

The  following  are  the  model’s  primary  simplifying  assumptions: 

•  Flat  Earth,  Flat  Terrain:  At  ranges  typical  for  SAL  guidance  (<25  km),  the  curvature  of  the 
earth  does  not  significantly  affect  the  model  geometry.  Flat  terrain  greatly  simplified  the 
model  development,  although  it  is  listed  as  a  potential  area  for  further  development  in 
section  7  of  this  report  (Path  Forward). 

•  Linear  Optics:  Because  semi-active  designator  lasers  are  typically  of  relatively  low  power, 
nonlinear  optical  effects  (e.g.,  thennal  blooming)  are  not  included  in  the  model. 

•  Simplified  Atmospheric  Transmission:  Atmospheric  attenuation  through  absorption  and 
reflection  is  incorporated  only  as  a  gross  power  loss  to  the  transmission.  It  is  assumed  that 
the  power  reflected  or  reemitted  from  the  transmission  path  to  the  seeker  is  negligible  (i.e. 
phenomena  such  as  backscatter  are  not  modeled). 

•  Lambertian  Reflection:  Lambertian  reflection  assumes  perfectly  diffuse  (matte)  reflective 
surfaces.  This  is  another  suggested  area  for  further  development  discussed  in  section  7. 

The  C++  class  implementation  of  the  model,  sSalSeeker,  can  be  executed  stand-alone  or  be 
embedded  into  a  larger  guided  projectile  model.  The  model  was  validated  using  the  Night 
Vision  Laser  Designator  (NVLaserD)  model  written  by  the  Night  Vision  and  Electronic  Sensors 
Directorate  (NVESD)  of  the  Communications-Electronics  Research,  Development,  and 
Engineering  Center  (CERDEC). 

The  report  is  organized  into  the  following  sections:  Section  2  describes  the  model  coordinate 
systems  and  major  entities  (bodies).  Section  3  explains  the  underlying  physics  of  the  laser 
transmission  algorithm.  Section  4  describes  the  seeker  guidance  algorithm.  Section  5  discusses 
the  model’s  C++  class  implementation.  Section  6  gives  the  results  of  the  model’s  validation,  and 
we  conclude  in  section  7  with  some  recommendations  for  the  use  and  further  development  of  the 
model. 
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2.  Model  Setup 


2.1  Global  Coordinate  System 

The  global  coordinate  system  is  Cartesian  and  fixed  to  the  earth’s  surface,  which  is  assumed  to 
be  an  infinite  plane.  The  global  coordinate  system  is  defined  by  three  orthogonal  unit  vectors: 
x,  y,  and  z.  The  x-y  plane  forms  the  earth’s  surface  (ground  plane),  where  z  =  0.  As  z  increases, 
the  height  above  the  earth’s  surface  increases.  Figure  1  showed  the  model  set-up  in  global 
coordinates. 

The  orientation  of  the  x  and  y  axes  and  the  coordinate  system  origin  are  defined  by  the  user  with 
the  following  conditions: 

1 .  The  x  and  y  axes  may  be  oriented  in  any  direction  on  the  ground  plane,  provided  they 
follow  the  right-hand  rule  (x  x  y  =  z). 

2.  The  origin  of  the  global  coordinate  system  may  be  situated  anywhere  on  the  ground  plane. 

The  global  coordinate  system  is  the  default  system  used  throughout  the  model.  If  not  specified, 
one  should  assume  these  coordinates.  A  major  purpose  of  the  global  coordinate  system  is  to 
define  the  relative  positions  and  orientations  of  the  model’s  bodies,  which  is  described  in  the 
next  section. 

2.2  Bodies  (Designator,  Target,  Seeker) 

The  bodies  involved  in  SAL  guidance  are  the  designator,  target,  and  seeker.  The  user  is 
responsible  to  input  all  body  positions  and  orientations,  with  the  exception  of  designator 
orientation.  From  the  orientations,  the  model  calculates  normalized  pointing  vectors  (unit 
vectors)  for  each  body. 

2.2.1  Designator 

The  designator  is  assumed  to  be  a  point-mass,  whose  coordinates  are  inputs  to  the  model.  The 
designator  pointing  vector,  d,  points  from  the  designator  location  toward  the  target’s  geometric 
center. 

2.2.2  Target 

The  target  is  a  rectangular  parallelepiped,  with  dimensions  of  length,  width,  and  height  (figure 
2).  The  target  is  assumed  to  always  stay  upright,  meaning  its  top  surface  stays  parallel  to  the 
ground  plane.  The  target  pointing  vector,  t,  points  from  the  target’s  geometric  center  toward  the 
center  of  the  target’s  front  surface  (shown  in  blue  in  the  figure).  The  target’s  right  and  left  sides 


* 

All  unit  vectors  will  be  denoted  by  the  A  symbol. 
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are  defined  by  the  respective  right-  and  left-handed  directions  when  facing  from  the  target’s 
geometric  center  towards  the  front  of  the  target. 

The  orientation  of  the  target  is  specified  by  the  user  through  a  target  rotation  angle,  x,  which  is 
the  angle  between  the  global  y  axis  and  the  projection  of  the  vector  pointing  away  from  the 
target’s  front  surface  onto  the  x-y  plane  measured  clockwise  as  seen  from  above  ( Bounds :  0  to 
2n). 

2.2.3  Seeker 

The  model’s  SAL  seeker  is  modeled  as  a  four-quadrant  IR  detector  mounted  in  the  nose  of  a 
projectile.  It  is  modeled  as  a  zero-thickness  circular  disc.  The  seeker  pointing  vector,  s,  aligns 
with  the  long  axis  of  the  projectile  and  points  from  the  center  of  the  seeker  normal  to  the  seeker’s 
plane  and  away  from  the  projectile.  The  seeker’s  FOV  is  the  angular  measure  of  the  cone 
centered  about  s  in  which  the  seeker  receives  radiation. 

The  user  specifies  s  through  azimuth  and  elevation  angles: 

•  Azimuth  (4>)  -  the  angle  between  the  global  x  axis  and  the  projection  of  s  onto  the  x-y 
plane.  Azimuth  is  positive  in  the  counter-clockwise  direction  as  seen  from  above.  Bounds: 
-7 r  to  K. 

•  Elevation  (9)  -  the  angle  between  the  global  x-y  plane  and  s.  Elevation  is  positive  when  s 
points  above  the  x-y  plane.  Bounds:  -n/2  to  n/2. 

The  seeker  is  divided  by  yaw  and  pitch  axes,  which  will  be  discussed  in  the  next  section. 
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2.3  Body-Fixed  Coordinate  System 


A  Cartesian  body-fixed  coordinate  system  (x',  y',  z')  is  used  to  define  vectors  relative  to  a  body 
pointing  vector,  b  (section  2.2).  In  this  -prime”  coordinate  system  used  by  Yager  (7),  the  origin 
is  set  at  the  geometric  center  of  the  local  body  (equation  1).  The  x'  axis  is  aligned  with  the  body 
pointing  vector,  b.  The  y'  axis  points  orthogonal  to  the  left  of  b  and  parallel  to  the  global  x-y  plane. 
The  z!  axis  is  fonned  orthogonal  to  the  x'  and  y'  axes  according  to  the  right-hand  rule  (figure  3). 


In  equation  fonn: 


x'  =  b. 

(1) 

'  —  z  xb. 

(2) 

-  x'  x  y' . 

(3) 

»  ‘ 
X 

\ 

r 

PS  Global  Coordinates  \ 

V 

n 

|  Body-Fixed  Coordinates  \ 

/  < 

f)  Body  Pointing  Vector  in  Global  Coordinates 

Defines  x' axis  in  Body-Fixed  Coordinates  \ 

V  Normalized  Vector  in  Body-Fixed  Coordinates 

To  be  rotated  into  Global  Coordinates  \ 

V  Normalized  Vector  in  Global  Coordinates  / 

X 

Figure  3.  Body-fixed  coordinates. 

In  the  prime  system,  a  vector,  v' ,  is  defined.  The  methodologies  to  define  v'  are  discussed  in 
sections  2.3. 1-2. 3. 3  of  this  report.  To  transform  v'  from  the  body-fixed  system  to  the  global 
coordinate  system,  we  translate  to  the  location  of  the  body  in  the  global  system  and  then  rotate  to 
align  the  x'  axis  with  the  body’s  pointing  vector,  b.  To  accomplish  this  rotation,  the  model  uses 
a  rotation  matrix  that  is  derived  by  Yager  ( 1 ).*  To  find  the  vector  in  global  coordinates  (fi),  the 
rotation  matrix,  A,  is  multiplied  by  v'\ 


* 

The  rotation  matrix  is  shown  in  a  simplified  form  used  for  unit  vectors. 
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v  —  Av' , 


(4) 


where 


A  = 


bx 

by 

bz 


-yby  -ybxbz 
-ybx  —ybybz 

0  - 

Y 


(5) 


where  bx  ,by  ,  and  bz  are  the  components  of  the  body  pointing  vector,  b,  in  global  coordinates 


and: 


i 


(6) 


2.3.1  Designator  Application 

For  the  designator,  body-fixed  coordinates  are  used  for  two  purposes: 

1 .  Adjusting  the  designator  pointing  vector  (d)  for  aim  error.  This  new  pointing  vector  is 
called  da[m. 

2.  Modeling  the  beam’s  solid  divergence  cone  that  is  centered  on  da im.  This  is  done 
stochastically  by  dividing  the  beam  into  a  set  of  rays.  Body- fixed  coordinates  are  used  to 
calculate  each  individual  ray  vector,  find. 

2.3. 1 . 1  Aim  Error.  To  model  aim  error,  d  is  perturbed  through  horizontal  and  vertical  perturbation 
angles,  ay  and  az\ 

Referring  to  figure  4,  the  perturbed  vector  is  expressed  in  the  body-fixed  coordinate  system  by 

adding  the  perturbations  to  the  original  pointing  vector. 

,  r 

z  Components  of  Modified  Vector  (  v  ) 


Figure  4.  Modified  vector  in  body-fixed  coordinates. 
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v'  =  1  ■  x'  +  tan(ay)  •  y'  +  tan(ay )  •  z' . 
To  normalize  the  perturbed  vector,  v' ,  it  is  divided  by  its  magnitude: 


(V) 


\v'\  —  Jl~^Xtan(ay) y^+~(tan(a^)y^. 


(8) 


(9) 


Finally,  we  transfonn  v'  into  the  global  coordinate  system  using  equation  4.  This  results  in  our 
new  designator  pointing  vector,  daim. 

2.3. 1.2  Solid  Beam  Divergence  Cone.  To  calculate  the  direction  of  each  ray  (find)  in  the  beam’s 
solid  divergence  cone,  we  define  a  perturbed  vector,  v' ,  in  the  same  way  for  aim  error  (equations 
7-9,  figure  5).  The  beam  is  assumed  to  be  circular,  so  that  ay  and  ay  vary  over  the  same  range. 
The  selection  of  ay  and  ay  is  done  through  random  draws  from  a  normal  distribution,  which  is 
discussed  in  section  3.1.1. 


Figure  5.  Calculation  of  individual  ray  vectors. 


2.3.2  Target  Application 

Because  of  the  target’s  simplified  geometry  and  rotations,  the  model  does  not  currently  use 
body-fixed  coordinates  for  the  target.  However,  if  the  model  were  updated  to  include  complex 
target  geometry  or  rotations,  body-fixed  coordinates  could  greatly  simplify  the  characterization 
of  target  surfaces. 

2.3.3  Seeker  Application 

For  the  seeker,  body-fixed  coordinates  are  used  to  find  the  yaw  and  pitch  axes  in  the  seeker 
plane.  In  the  body-fixed  coordinate  system,  the  pitch  axis  is  the  y'  axis,  and  the  yaw  axis  is  the  z' 
axis  (figure  6). 
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Thus,  to  find  the  pitch  and  yaw  axes  in  global  coordinates,  we  first  set  the  vector  v'  equal  to  each 
axis  in  body-fixed  coordinates: 

vpitcn  =  y'  =  (0,1,0).  (10) 

v'  =  z'  =  (0,0,1).  (11) 

After  rotating  back  to  global  coordinates  using  the  rotation  matrix  (A),  we  have  the  pitch  and 
yaw  axes  in  global  coordinates. 


3.  Laser  Transmission  Model 


The  laser  transmission  algorithm  follows  four  successive  stages  (figure  7).  In  each  stage,  the 
model  calculates  the  beam’s  transmission  path  and  power  loss: 

•  Stage  1:  Atmospheric  Transmission  (Designator  to  Target) 

•  Stage  2:  Target  Reflection 

•  Stage  3:  Atmospheric  Transmission  (Target  to  Seeker)  7 

l  Combined  into  one 

•  Stage  4:  Seeker  Reception  J  stage  in  this  report 
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Figure  7.  Laser  transmission  stages. 


3.1  Stage  1:  Atmospheric  Transmission 

3.1.1  Beam  Divergence 

Prior  to  modeling  the  beam  divergence,  the  beam’s  centerline  must  first  be  defined.  The 
designator  is  assumed  to  always  aim  for  the  geometric  center  of  the  target,  and  this  pointing 
vector  is  perturbed  in  random  horizontal  and  vertical  directions  to  account  for  designator  error, 
as  described  in  section  2.3.1 .  The  perturbation  angles  are  drawn  from  normal  distributions  with 
standard  deviations  aUyi  and  oa  , .  This  results  in  the  vector  daim. 

As  the  beam  emerges  from  the  designator,  it  diverges  along  its  transmission  path.  Assuming  the 
beam  to  have  a  Gaussian  profile  across  its  transverse  axis,  the  divergence,  5,  is  commonly 
defined  to  be  the  half  angle  corresponding  to  the  location  along  the  transverse  axis  where  the 
intensity*  drops  to  1/e2  times  the  intensity  at  the  beam  centerline  (figure  8)  (2).  At  any  range 
from  the  designator,  the  diameter  of  the  cone  swept  out  at  this  point  on  the  transverse  axis  is 


Intensity  is  power  per  unit  solid  angle,  and  it  will  be  discussed  in  section  3.1.2. 
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referred  to  as  the  beam  diameter.*  This  occurs  at  two  standard  deviations  from  the  beam 
centerline  intensity,  meaning  that  -95%  of  the  laser  beam’s  total  power  is  within  the  cone  of  the 
given  divergence  angle. 


Figure  8.  Gaussian  beam  with  1/e2  beam  diameter. 


Next,  the  beam  is  divided  stochastically  into  a  set  of  rays,  with  each  ray  carrying  its  individual 
portion  of  the  total  power.  Each  ray’s  direction  (find)  is  determined  using  body-fixed 
coordinates,  as  described  in  section  2.3.1.  The  perturbation  angles,  ay  and  az',  are  drawn  from  a 
normal  distribution  with  a  standard  deviation  based  on  the  beam  divergence  angle  (cra  ,  ,  =  as). 

Recalling  that  the  beam  diameter  corresponding  to  8  is  defined  by  2atransverse5  which  varies 
closely  to  2 os,  we  solve  for  er5f 


8  =  2  os. 

(12) 

s 

(13) 

3.1.2  Attenuation 


The  first  part  of  atmospheric  transmission  was  determining  the  direction  of  the  beam,  which  was 
accomplished  through  a  division  into  a  set  of  rays,  each  with  its  own  direction,  find.  The  second 
part  of  atmospheric  transmission  is  determining  the  power  loss  for  each  ray  on  its  path  to  the 
target.  The  starting  power  for  each  ray  (())ray)  is  determined  by  dividing  the  beam’s  pulsed  power 
by  the  total  number  of  rays  (n),  and  multiplying  by  a  designator  efficiency  coefficient 


(^efficiency)- 


* 

There  are  alternative  methods  for  defining  beam  divergence.  It  is  sometimes  defined  as  the  full  cone  angle  of  the  beam, 
instead  of  the  half  angle  that  we  are  assuming.  In  addition,  the  beam  width  is  sometimes  defined  according  to  the  full-width  at 
half  maximum  method  (FWHM). 

^This  is  a  small  angle  approximation,  as  the  divergence  by  definition  produces  a  nonnal  distribution  in  distance  across  the 
beam’s  transverse  axis  (CTtransverse),  and  not  in  the  beam’s  divergence  angle.  The  tangent  of  the  divergence  angles,  ay  and  az,  is 
linked  to  the  transverse  distance.  For  small  angles,  tana  =  a,  and  the  divergence  angles  approximately  follow  a  normal 
distribution. 
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efficiency- 


(14) 


4*ray 


^pulse 


•  d 


Atmospheric  attenuation  is  the  exponential  decrease  in  beam  intensity  as  it  transmits  through  the 
atmosphere.  The  Beer-Lambert  Law  characterizes  this  attenuation: 

l  =  l0e~kx,  (15) 


where: 


I  =  Attenuated  Intensity  at  distance  x, 

Io  =  Initial  Intensity, 
k  =  Attenuation  Coefficient,  and 
x  =  Path  Length. 

Beam  intensity,  also  known  as  radiant  intensity,  is  defined  as  the  power  (tj))  per  unit  solid  angle 
(Q)  subtended  by  the  beam.  The  solid  angle  is  the  surface  area  subtended  on  a  sphere  of  radius  r 
(figure  9).  It  is  measured  in  steradians  (sr),  where  the  sphere  represents  4n  sr. 


Figure  9.  Solid  angle. 


Thus,  the  solid  angle  is  calculated  by  dividing  the  subtended  surface  area  (S)  by  the  square  of  the 
sphere’s  radius  (r): 

«=£■  (1«) 

In  the  case  of  a  laser  beam,  r  is  the  range  traversed  by  the  beam  measured  along  its  centerline. 
Because  the  model  does  not  consider  nonlinear  optical  phenomena,  the  beam  divergence  remains 
constant  throughout  the  transmission.  This  means  that  the  solid  angle  also  remains  constant,  and 
therefore,  in  this  case,  we  can  generalize  the  Beer-Lambert  Law  for  power  (see  figure  10  and 
equations  17-20). 
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For  small  8, 


Thus, 


nx 


Si  =  it  •  (rx  •  5)2,  S2  s  rt  •  (r2  5)2. 


Si  _  Ti-(rr5)2 

ri2  ri2 


-rt  •  52,  n2 


S2  Ti-(r2-5)2  _ 

r22  r22 


Hi  =  n2 . 


(17) 

(18) 


(19) 


Therefore,  for  a  beam  of  constant  divergence,  we  can  generalize  the  Beer-Lambert  Law  for  beam 
power,  (j): 

t  _  ^  _  jfo  -kx 

n  n0  ’ 

$  =  <$>oe~kx.  (20) 


Because  of  the  complexity  and  variability  of  the  earth’s  atmosphere,  the  attenuation  coefficient 
(k)  varies  according  to  many  factors,  and  it  needs  to  be  calculated  for  each  scenario.  Equation  21 
approximates  the  attenuation  coefficient  by  summing  four  major  components: 


k  k]y[0]ecu]ar  Absorption  T  kMoleeiilar  Scattering  “E  ^Aerosol  Absorption  T  k^erosoi  Seatteriag- 


(21) 


Molecular,  in  the  context  of  equation  21,  refers  to  atmospheric  particles  larger  than  electrons  but 
smaller  than  X,  the  laser  wavelength.  Similarly,  aerosol  refers  to  particles  that  have  a  size 
comparable  to  X  (2).  Given  these  four  attenuating  components,  the  model  utilizes  lookup  tables 
to  determine  the  attenuation  coefficient  (3).  The  tables  break  up  the  atmosphere  into  one- 
kilometer-deep  altitude  segments.  The  look-up  tables  require  the  following  infonnation: 
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•  Start  Height,  End  Height,  Path  Length  -  if  the  path  traverses  multiple  altitudes,  the 

model  breaks  up  the  path  and  steps  through  each  altitude  segment 

•  Laser  Type  -  characterized  by  the  wavelength.  The  look-up  tables  currently  handle  two 

laser  types: 

Nd:YAG  laser  (1.06*  pm)  -  the  most  popular  wavelength  currently  for  range-finding 
and  designation,  non-eye-safe. 

Er:Glass  laser  (1.54  pm)  -  an  eye-safe  wavelength  of  interest. 

•  Visibility  -  affecting  aerosol  absorption  and  scattering 

Clear:  23  km 
Hazy:  5  km 

•  Latitude  -  affecting  molecular  absorption  and  scattering 

Tropics:  0°  to  23.5°  and  0°  to -23.5° 

Mid-Latitudes:  23.5°  to  50°  and  -23.5°  to  -50° 

Sub-Arctic:  50°  to  70°  and  -50°  to  -70° 

•  Season  -  affecting  molecular  absorption  and  scattering 

Summer:  March  22  -  September  2 1 
Winter:  September  22  -  March  2 1 

Thus,  using  equations  20  and  21,  the  model  determines  the  attenuated  power  of  the  ray  where  it 
intersects  the  target  surface  (  c|)target).  The  distance  from  the  designator  to  surface  intersection 
point  in  equation  20  (x),  is  determined  by  ray  tracing  in  section  3.2.1. 


3.2  Stage  2:  Target  Reflection 


3.2.1  Ray  Projection  onto  Target 


The  first  element  of  target  reflection  is  determining  if  and  where  each  ray  hits  the  target.  To  this 
end,  we  employ  ray  tracing,  which  determines  the  first  plane  of  interest  that  is  intersected  by  the 
ray  and  the  intersection  point  on  that  plane.  Wikipedia  describes  the  intersection  in  matrix 
notation  ( 4 ): 


-t- 

- -y*  v  -v/*  y  y*  _  y 

xa  -A-fo  A-i  Aq  A2  Xq 

-1 

■xa  -  x0- 

u 

— 

ya  -  jb  yi-  y  o  y-i  -  y0 

ya  -  y0 

(22) 

-V- 

za  ~  zb  zl~  z0  z2  ~  z0. 

za  z0. 

NctGlass,  a  laser  of  nearly  identical  wavelength  to  Nd:YAG,  was  substituted  in  the  attenuation  lookup  tables. 
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where 


P0  (x0,  yo,  z0)  =  reference  point  on  plane, 

Pi  (xlf  yi,  z\ )  =  second  point  on  plane  (defining  first  direction  in  plane  relative  to  PG), 

P2  (x2,  y2,  z2)  =  third  point  on  plane  (defining  second  direction  in  plane  relative  to  P0), 

la  (xa,  ya,  za)  =  starting  point  of  the  ray, 

lb  (xb,  yb,  Zb)  =  second  point  on  ray,  defining  the  direction, 

t  =  distance  between  ray  start  point  and  the  plane, 

u  =  distance  in  plane  from  Po  to  ray  intersection  in  first  direction,  and 

v  =  distance  in  plane  from  Po  to  ray  intersection  in  second  direction. 

In  the  model,  these  elements  are  defined  in  global  coordinates  (figure  11): 

Input 

P0  =  center  point  of  surface. 

Pi  =  center  point  of  left  edge  of  the  surface. 

P2  =  center  point  of  top  edge  of  the  surface. 
la  =  designator  location. 

lb  -  la  =  individual  ray  vector  from  designator  ( find ). 

Output 

t  =  distance  between  designator  and  surface  intersected, 
u  =  distance  from  center  of  surface  to  ray  intersection  in  u  direction, 
v  =  distance  from  center  of  surface  to  ray  intersection  in  v  direction, 
lb  =  ray  hit  point  on  plane  =  P0  +  uu  +  w. 
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To  see  if  the  ray  intersects  a  finite  surface,  u  and  v  are  compared  to  the  dimensions  of  the 
surface.  For  the  parallelepiped  target,  this  is  half  of  the  target  length,  width,  or  height,  depending 
on  the  surface.  If  u  and  v  are  both  less  than  one  half  of  the  dimension,  there  exists  an 
intersection,  unless  the  ray  intersects  another  surface  first.  Thus,  the  model  implements  the  ray¬ 
tracing  routine  for  all  seven  surfaces,  checking  to  see  what  surface(s)  the  plane  intersects,  and,  of 
those,  choosing  the  one  with  the  shortest  ray  length,  t.  The  first  five  surfaces  tested  form  a 
rectangular  parallelepiped  target:  Left  Side,  Front,  Right  Side,  Back,  and  Top.*  The  sixth 
surface  is  the  ground  plane.  Ray  intersection  with  the  ground  plane  indicates  underspill  or 
overspill.  The  final  surface  that  the  ray-tracing  routine  checks  is  the  seeker  itself.  This  tests  the 
rare  situation  in  which  the  laser  ray  is  pointing  directly  into  the  seeker  from  the  designator. 


* 

The  target’s  bottom  surface  is  not  tested  in  the  ray  tracing  routine,  although  it  could  be  added  to  the  model  in  the  future. 
This  capability  could  be  useful  if  the  target  was  in  the  air,  and  the  projectile/seeker  had  the  capability  to  fly  upwards  towards  the 
target. 
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3.2.2  Surface  Reflection 


Once  the  ray  hits  a  surface,  it  will  either  be  absorbed  or  reflected.  To  model  the  power  lost  at  the 
surface  due  to  absorption,  the  model  employs  a  simple  -target  reflectivity”  multiplier  (Reflectivity), 
which  can  be  varied  for  each  surface: 

^target, reflected  —  4*target, received  '  ^reflectivity.  (23) 

Several  studies  have  investigated  the  reflective  properties  of  different  targets,  which  are  a  result 
of  both  surface  composition  and  laser  wavelength  (5). 

There  are  two  primary  types  of  reflection  (figure  12): 

1 .  Specular  Reflection  -  mirror-like  reflection,  where  the  angle  of  incidence  (9;) 
equals  the  angle  of  reflection  (0r). 

2.  Lambertian  Reflection  -  diffuse  reflection,  where  the  reflected  energy  is  scattered 
in  all  directions  regardless  of  angle  of  incidence. 


Figure  12.  Specular  vs.  Lambertian  reflection. 

Real-world  scenarios  generally  have  elements  of  both  types  of  reflection,  but  Lambertian 
reflection  is  closer  to  actual  behavior,  and  it  is  used  exclusively  in  the  model.  Thus,  the  model 
assumes  that  the  energy  reflects  off  of  the  target  in  all  directions  in  a  180°  hemisphere  normal  to 
the  surface  of  reflection.  This  will  be  described  later  in  section  3.3.2. 

3.3  Stages  3  and  4:  Atmospheric  Transmission  and  Seeker  Reception 

3.3.1  Ray  Projection  Into  Seeker 

After  target  reflection,  the  model  has  calculated  the  following  data  for  each  ray: 

•  Coordinates  of  ray  surface  intersection  (lb) 

•  Surface  the  ray  intersects  (Left  Side,  Front,  Right  Side,  Back,  Top,  Ground,*  Direct-to- 
Seeker,  or  no  surface  hit) 


If  the  ray  hits  the  ground,  the  model  tests  whether  the  ray  is  obscured  from  the  seeker  by  the  target. 


16 


Ray  Power  ((jltarget, reflected) 


Using  these  data,  we  can  now  find  if  and  where  the  rays  project  into  the  seeker’s  four-quadrant, 
IR  sensor.  To  find  if  a  ray  projects  into  the  seeker,  it  must  pass  two  tests  (refer  to  figure  1): 

1 .  Correction  Angle  (p)  between  seeker  heading  (s)  and  the  vector  from  the  seeker  to  the  ray 
hit  point  ( st )  must  be  less  than  or  equal  to  the  seeker  FOV  (Ps). 

2.  Off- Angle  (a)  between  the  surface  normal  vector  of  ray  hit  (n)  and  the  vector  from  the  ray 
hit  point  to  the  seeker  (ts)  must  be  less  than  90°. 


These  angles  are  detennined  using  the  dot  product  of  the  two  vectors: 


*  = cos_1  (jfrii)  -  P*  ■  <24> 


If  the  ray  does  not  pass  both  of  these  tests,  it  does  not  project  into  the  seeker’s  sensor. 


If  the  ray  passes  both  tests,  the  model  determines  which  of  the  four  seeker  quadrants  the  ray 
projects  into.  To  do  so,  the  model  again  uses  ray  tracing.  In  this  instance,  the  ray  originating  at 
the  target  intersection  (h)  is  pointed  in  the  negative  direction  of  s,  and  it  is  determined  where  it 
intersects  the  seeker  plane  (figure  13).  The  pitch  and  yaw  axes  forming  the  seeker  plane  are 


r  ti 


calculated  using  body-fixed  coordinates  (section  2.3.3).  The  ray  tracing  produces  a 


it 

-v- 


output 


vector  that  locates  the  ray’s  intersection  point  with  the  seeker: 


t  =  ray  distance  to  seeker  plane. 


u  =  distance  from  seeker  center  to  ray  intersection  in  the  yaw  direction. 


v  =  distance  from  seeker  center  to  ray  intersection  in  the  pitch  direction. 


A  seeker  intersection  in  the  positive  yaw  direction  (+u)  indicates  that  the  seeker  must  rotate 
positively  about  the  yaw  axis  (to  the  left  in  the  FOV)  to  point  at  (refer  to  figure  6).  The 
positive  pitch  direction  (+v)  is  similarly  linked  to  a  positive  rotation  about  the  pitch  axis  (to  the 
top  in  the  FOV).  Because  we  already  know  that  the  ray  projects  into  the  seeker’s  sensor,  we  only 
need  to  know  whether  u  and  v  are  positive  or  negative  to  determine  what  quadrant  the  ray 
projects  into. 
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Figure  13.  Ray  intersection  with  seeker  plane. 

3.3.2  Power  Received  by  Seeker 

After  recording  what  quadrant  the  ray  projects  into,  the  model  determines  the  ray’s  individual 
power  contribution  to  the  quadrant.  To  determine  the  amount  of  power  received  by  the  seeker 
for  each  ray,  the  model  uses  Lambert’s  Cosine  Law.  The  law  states  that  the  radiant  intensity  I 
(power  per  steradian)  received  from  a  perfectly  Lambertian  surface  is  proportional  to  the  cosine 
of  the  angle,  a,  between  the  observer’s  line  of  sight  and  the  surface  normal  (figure  1,  equation  25). 

In  reference  6,  McCartney  describes  how  to  calculate  the  radiant  flux  0  (power)  across  a 
hemisphere  by  integrating  the  radiant  intensity,  I  (figure  14).  In  the  hemisphere,  a  represents  the 
90°  complement  to  elevation,  and  to  represents  azimuth.*  From  the  conservation  of  energy,  we 
know  that  the  total  power  in  the  system  passing  through  the  hemisphere  is  equal  to  the  total 
power  reflecting  off  of  the  target  (neglecting  attenuation  losses,  which  will  be  factored  in  later). 
Therefore,  if  we  integrate  over  the  hemisphere  of  intensity  Ipeak'COsa,  we  get  the  total  power 
passing  through  the  hemisphere,  0  (equations  26-27).  Because  we  know  0  as  the  reflected 
power  off  of  the  target,  we  can  solve  for  the  hemisphere’s  peak  radiant  intensity,  Ipeak  (equation 
28).  Finally,  using  Lambert’s  Cosine  Law,  we  determine  the  radiant  intensity  at  the  seeker’s 
position  on  the  hemisphere  (equation  29). 


* 

Note:  a  and  co  are  different  from  the  seeker  azimuth  (<)))  and  elevation  (0).  a  and  co  characterize  the  hemisphere  extending 
from  the  designated  surface  and  the  seeker’s  position  on  this  hemisphere.  (])  and  0  describe  the  orientation  of  the  seeker  in  the 
global  coordinate  system. 
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Figure  14.  Integration  of  radiant  intensity  over  a  hemisphere. 


0  =  -C^2  Ipeakcosa(sina)d(jOda. 

Substituting  a  double  angle: 

r-n/2  r2n  l 

0  =  Ipeak  -sin(2a)duda 

J  o  Jo  z 

rn/2 

<p  —  n  •  Ipeak  I  sin(2a)da 
Jo 

0  —  Ft  ‘  Ipeak- 

Therefore,  solving  for  the  peak  radiant  intensity  in  the  hemisphere: 

T  _  ^target; ref lected 

Ipeak  ~  ~  • 


(26) 


(27) 


(28) 


Using  Lambert’s  Cosine  Law,  we  determine  the  radiant  intensity  at  the  seeker’s  location  on  the 
hemisphere: 


/ 


I  seeker  ^pea/c  '  COS 

41  tar  get, re  fleeted 


seeker 


cos  a. 


(29) 


To  detennine  the  power  at  the  seeker  (before  factoring  in  other  losses),  we  multiply  Iseeker  by  the 
solid  angle  subtended  by  the  seeker  on  the  hemisphere.  In  the  model,  the  hemisphere  extends 
from  the  laser  target  intersection  point  nonnal  to  the  surface  of  intersection.  The  radius  of  the 
hemisphere,  r,  is  the  range  from  the  laser-target  intersection  point  to  the  geometric  center  of  the 
seeker.  To  approximate  the  number  of  steradians  the  seeker  subtends  on  the  hemisphere 
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2 

(assuming  rseeker «  r),  the  seeker  area,  Sseeker,  is  divided  by  r  and  multiplied  by  the  cosine  of  the 
seeker  correction  angle,  p  (figure  1,  equation  16): 

ft seeker  =  •  COSp  (30) 


Thus,  including  attenuation  from  target  hit  to  seeker  (LAtten),  and  seeker  efficiency  (seffiCiency),  the 
laser  power  scattered  from  the  target  into  the  seeker  per  ray  is 


tfiseeker  I seeker  ’  ft seeker  '  L Atten  '  ^efficiency 

i  Vtarget; ref  lected  '  cos  a  '  $ seeker  *  C0SM^ '  ^Atten  ‘sefficiency 

0 seeker  =  - ~2 - 


(31) 


Equation  3 1  is  calculated  for  each  ray.  The  model  sums  the  individual  ray  contributions  for  each 
quadrant  to  find  the  total  power  received  in  all  four  seeker  quadrants.  The  next  section  discusses 
how  the  model  uses  these  data  to  produce  guidance  signals. 


4.  Seeker  Guidance  Model 


In  the  last  section,  the  model  calculated  the  beam’s  transmission  path  and  power  loss,  resulting  in 
a  power  distribution  at  the  seeker.  In  this  section,  we  discuss  how  the  model  interprets  these  data 
to  produce  projectile  guidance. 

4.1  Target  Encounter  and  Detect 

There  exists  a  progressive  ladder  in  the  target  knowledge  of  a  terminal  seeker  (figure  15). 


Encounter  Detect  Classify  Identify 

•Target  is  within  *1  see  »l  can  put  it  into  I^HI  *1  know  it  is  a 

my  field  of  something  a  group,  such  legitimate 

view  (not  ^L^B  standing  out  ^L^B  as  vehicles,  ^L^B  target 

directly  from  the  buildings,  or  Kj&Sg 

sensed)  background  natural  &f|||j| 

Bftl  SB  features  MB 


Figure  15.  Target  knowledge  ladder. 

For  the  SAL  seeker,  the  knowledge  ladder  is  modeled  in  the  following  way: 

•  Encounter  is  determined  in  the  model  by  calculating  the  angle  the  seeker  needs  to  rotate 
through  to  point  at  the  target’s  center.  If  the  angle  is  less  than  or  equal  to  the  FOV,  it 
returns  true.  Encounter  is  internal  to  this  model,  as  it  is  a  geometric  calculation  and  is  not 
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sensed  by  a  real  projectile.  It  is  a  useful  calculation,  however,  to  determine  the  point  in 
the  ballistic  trajectory  at  which  the  seeker  is  pointing  close  enough  to  the  target  to 
encounter,  which  is  the  first  step  toward  detection. 

•  Detection  occurs  in  the  model  when  the  laser  power  the  seeker  receives  is  above  the  S/N 
threshold  necessary  to  distinguish  the  signal.  The  noise  is  determined  using  a 
background  noise  multiplier,  which  assumes  that  the  seeker  noise  increases  linearly  with 
seeker  area  and  solid  angle  subtended  by  the  FOV.  This  also  assumes  that  the  seeker 
internal  noise  is  insignificant  compared  to  the  external  background  noise,  as  internal 
noise  does  not  scale  with  seeker  area  or  solid  angle.  The  model’s  default  multiplier  value 
is  based  on  a  commercial  seeker  with  minimum  detectable  signal  irradiance  (power  per 
unit  seeker  area)  of  35nW/cm  ,  FOV  half  angle  of  4.5°,  and  an  assumed  S/N  ratio  of  7 
(7). 

•  Classify  and  Identify  are  combined  in  the  final  step,  and  occur  when  the  SAL  seeker 
examines  the  signal’s  pulse  width  and  pulse  frequency  to  weed  out  false  signals.  Because 
the  model  does  not  contain  a  progressive  time  element,  these  effects  are  not  modeled,  but 
could  easily  be  added  upon  integration  into  a  larger  guided  trajectory  program. 

4.2  Guidance  Updates 

After  the  seeker  detects,  it  needs  to  process  the  information  into  guidance  signals.*  It  does  this 
by  comparing  the  signals  in  each  of  the  four  sections  of  the  detector.  Figure  16  shows  the  four- 
quadrant  detector  used  in  the  model. 


Figure  16.  Four-quadrant  laser  detector. 


While  the  model  always  computes  guidance  signals  as  a  convenience  to  the  user,  they  should  not  be  considered  reliable 
unless  also  accompanied  by  a  positive  detect  calculation. 
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If  Pi  is  the  power  received  in  quadrant  1,  and  likewise  for  the  other  quadrants,  Hubbard  describes 
how  the  spots’  location  can  be  approximated  in  the  following  form  (5): 


Ay 

„  (Pi+p2)— (P3+P4) 

->  Pitch  Guidance, 

(32) 

r 

P  total 

Ax 

„  (Pl+P3)-(P2+P4) 

->  Yaw  Guidance, 

(33) 

r 

P  total 

where  Ptotal  =  P1 

+  P2  +  P3+  P4. 

(34) 

The  pitch  and  yaw  guidance  range  in  value  from  -1  to  +1,  which  mean  full  maneuvers  in  the 
negative  and  positive  directions,  respectively.  In  calculating  guidance,  the  model  assumes  a 
non-spun  projectile,  although  spin  could  be  accommodated  in  a  fairly  simple  manner  upon 
integration  into  a  larger,  time-dependent  model. 

The  model  also  returns  a  polar  representation  of  where  the  power  is  received  on  the  sensor: 

•  Angle  (0):  0  to  2n  rad,  0  =  -Yaw  direction,  angle  increases  in  the  counter-clockwise 
(CCW)  direction. 

•  Magnitude:  0  to  1 . 

The  angle  detennines  how  to  combine  pitch  and  yaw  commands  and  the  magnitude  detennines 
how  much  to  maneuver  in  the  prescribed  direction.  The  polar  output  is  currently  a  repackaged 
version  of  the  preceding  pitch  and  yaw  guidance,  but  could  be  a  better  representation  to  use  in 
the  future  for  imaging  sensors,  in  order  to  track  multiple  data  points  within  a  single  FOV. 


5.  C++  Implementation 


5.1  Input  Variables 

The  input  variables  (table  1)  usually  change  throughout  an  analysis,  and  are  defined  in  the 
main.cpp  file. 
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Table  1.  Input  variables. 


Locations 

Type 

Units 

Bounds 

Description 

d_geo 

struct  Point3D 
Utilities  class 
double 

m 

z  >  =  0 

Designator  Location:  Global  Coordinates, 

Point  Mass 

t_geo 

Point3D 

m 

z  >  =  0 

Target  Location:  Global  Coordinates, 

Geometric  Center 

Due  to  flat  terrain,  the  ground  plane  is  located  at 
z  =  0.  If  the  target  is  ground-based,  its  z 
component  must  be  given  accordingly  (z  =  0  + 

1/2  target  height). 

s_geo 

Point3D 

m 

z  >  =  0 

Seeker  Location 

Global  Coordinates,  Geometric  Center 

Orientations 

Type 

Units 

Bounds 

Description 

s  orientation 

struct  Point2D 
Utilities  class 
double 

rad 

Az:  -7i  to  7i 

El:  -jt/2  to  tt/2 

Seeker  Orientation;  Globed  Coordinates 
.X  (Azimuth,  <)>)  -  the  angle  between  the  x  axis 
and  the  projection  of  s  onto  the  x-y  plane. 

Azimuth  is  positive  in  the  counter-clockwise 
direction  as  seen  from  above. 

.Y  (Elevation,  0)  -  the  angle  between  the  x-y 
plane  and  s.  Elevation  is  positive  when  s 
pointins  above  the  x-y  plane. 

t  rotation 

double 

rad 

0-2  7i 

Target  Rotation;  Global  Coordinates 
t  =  the  angle  between  the  +y  axis  and  the 
projection  of  the  target  heading  vector,  t,  onto 
the  x-y  plane  measured  clockwise  as  seen  from 
above,  t  points  from  the  target  geometric  center 
towards  the  center  of  the  front  face  of  the  target. 

Statistics 

Type 

Units 

Bounds 

Description 

seed 

int 

— 

— 

Seed  for  random  number  generator 

5.2  Input  Parameters 

The  parameters  (table  2)  tune  the  analysis,  but  generally  do  not  vary  during  a  simulation.  They 
can  be  set  in  the  main.cpp  fde,  or  through  an  input  fde  (see  appendices  A  and  B,  main.cpp  #2). 
If  no  change  in  a  parameter  is  detected,  the  program  uses  the  default  parameter  value,  which  is 
automatically  set  using  the  internal  function  SetDefaultsQ  when  the  class  is  first  instantiated. 
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Table  2.  Input  parameters. 


Weather 

Type 

Units 

Bounds 

Default 

Description 

sal  weather 

struct 

weather 

int 

— 

Seas:  0-2 
Lat:  0-1 
Vis:  0-1 

Seas:  0 
Lat:  1 

Vis:  0 

.season 

0  =  summer  (March  22  -  September  21) 

1  =  winter  (September  22  -  March  21) 
.latitude 

0  =  tropics  (0°  to  23.5°  ±) 

1  =  mid-latitudes  (23.5°  to  50°  ±) 

2  =  sub-arctic  (50°  to  70°  ±) 

.visibility 

0  =  clear  (23  km) 

1  =  hazy  (5  km) 

Seeker 

Type 

Units 

Bounds 

Default 

Description 

s_fov 

double 

rad 

0  -  nil 

0.079  (4.5°) 

Seeker  field  of  view  (half  angle  of  cone) 

saperture 

diameter 

double 

m 

>0 

0.060 

Seeker  aperture  diameter 

s_background 

double 

WrrUfi'1 

>  =  0 

0.010 

Seeker  background  noise  multiplier 

s  signalnoise 
threshold 

double 

— 

>0 

7.0 

Ratio  of  SAL  signal  to  background  noise  to 
be  able  to  detect 

s  efficiency 

double 

— 

0-1 

0.95 

Seeker  loss  coefficient 

Designator 

Type 

Units 

Bounds 

Default 

Description 

d_lasertype 

int 

— 

0-1 

0 

0  =  1.06  pm  (Nd:Glass) 

1  =  1.536  pm  (Er:  Glass) 

d_raycount 

int 

— 

=  1 

10000 

No.  of  rays  to  divide  laser  pulse  into 

ddivergence 

double 

rad 

>  =  0 

3E-4 

Divergence  of  laser  beam  (half  angle) 

d_pulse_energy 

double 

J 

>0 

80E-3 

Starting  designator  energy/pulse 

d_pulse_frequency 

double 

Hz 

>0 

10 

No.  ofpulses/second 

d_pulse_duration 

double 

s 

>0 

15E-9 

Pulse  length 

d_efficiency 

double 

- 

0-1 

0.95 

Designator  loss  coefficient 

d  h  error 

double 

rad 

>  =  0 

IE-4 

Designator  pointing  error  (horizontal,  SD) 

d_v_error 

double 

rad 

>  =  0 

IE-4 

Designator  pointing  error  (vertical,  SD) 

Target 

Type 

Units 

Bounds 

Default 

Description 

t  size 

struct 

Point3D 

Utilities 

class 

double 

m 

>0 

6.4,2.3,23 

Target  Dimensions 
.X  =  length 
.Y  =  width 
.Z  =  height 

t  reflect 

array[int] 

— 

0-1 

0.4  for  all 

Target  Reflectivity 

[0]  =  left  side,  [1]  =  front,  [2]  =  right  side 
[3]  =  back,  [4]  =  top,  [5]  =  ground  plane 
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5.3  Output 

The  user  only  needs  to  call  one  function,  S  SALSEEKER,  which  performs  all  necessary 
calculations  and  outputs  all  results  to  a  class-defined  struct,  — skmtputs”  (see  table  3). 


Table  3.  The  S  SALSEEKER  output  function. 


Type 

Units 

Description 

SSALSEEKER 
(d_geo,  t_geo,  s_geo, 
s  orientation, 
trotation,  seed) 

struct  saloutputs 

— 

Main  calculation  function 

OUTPUT 

.encounter  (bool):  target  is  within  seeker  field  of  view 
.detect  (bool):  S/N  above  threshold 

.sal_signals  (struct):  magnitude  of  peak  laser  power  sensed 
in  each  quadrant  (W) 

.ql  =  Positive  Pitch,  Positive  Yaw 
.q2  =  Positive  Pitch,  Negative  Yaw 
.q3  =  Negative  Pitch,  Positive  Yaw 
.q4  =  Negative  Pitch,  Negative  Yaw 
.actuator_signals  (struct):  lifting  and  turning  guidance 
based  on  sal_signals  quadrant  values 
.Pitch,  .Yaw  (-1  to  1) 

-1  is  full  maneuver  in  negative  direction  (-yaw,  -pitch) 

+1  is  full  maneuver  in  the  positive  direction. 

Partial  maneuvers  for  numbers  in  between  -1  and  +1. 
.Theta,  .Mag  (polar  representation  of  signal  in  FOV) 
.Theta  =  Polar  direction  of  signal  center  (0-271  rad) 

0  rad  =  -  yaw  direction,  rotate  CCW) 

.Mag  -  Polar  magnitude  of  direction  vector  (0-1) 

6.  Validation 


All  validation  tests  were  run  using  the  default  parameters  unless  otherwise  noted  (default 
parameters  listed  in  section  5.2). 

6.1  Power  Loss 

The  NVLaserD  model  was  used  to  validate  the  power  drop  of  the  laser  across  a  given  distance 
(9)  (see  table  4).  The  beam  divergence  (8)  was  varied  from  0  to  2.4  milliradians  (tnrad)  and  for 
each  divergence,  the  seeker  angle  from  the  target  surface  was  tested  at  0°  and  45°.  When  8 
equaled  0,  the  results  agreed  within  1%.  As  8  increased,  the  results  originally  differed.  This 
disagreement  is  most  likely  explained  by  a  difference  in  the  way  8  is  defined  between  models. 
Recall  equation  13  ( as  =  8/2),  which  assumed  a  8  corresponding  to  a  beam  radius  of  2  a  from 
the  peak  intensity.  If  8  was  instead  assumed  to  correspond  to  a  1  -ct  beam  radius,  the  resulting 
beam  spread  would  be  twice  as  great  as  the  2-a  beam  spread  (figure  17). 
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Table  4.  Energy  drop  comparison  between  sSalSeeker  and  NVLaserD. 


Test  Inputs 

TestI 

Test  2 

Test  3 

Test  4 

Designator  Pulsed  Energy  (mj) 

80 

80 

80 

80 

80 

80 

80 

80 

Designator  Pulse  Width  (ns) 

15 

15 

15 

15 

15 

15 

15 

15 

Designator  Pulsed  Power  (Peak  -  MW) 

5.3 

5.3 

5.3 

5.3 

5.3 

5.3 

5.3 

5.3 

Designator  Error  (mrad) 

0 

0 

0 

0 

0 

0 

0 

0 

Target  Rotation  (°) 

0 

0 

0 

0 

0 

0 

0 

0 

Range  (Designator  to  Target,  km) 

2 

2 

2 

2 

2 

2 

2 

2 

Range  (Seeker  to  Target,  km) 

5 

5 

5 

5 

5 

5 

5 

5 

k  (attenuation  coefficient) 

0.0581 

0.0581 

0.0581 

0.0581 

0.0581 

0.0581 

0.0581 

0.0581 

5  (divergence  -  mrad)  -  s SalSeeker 

0 

0 

0.6 

0.6 

1.2 

1.2 

2.4 

2.4 

8  (divergence  -  mrad)  -  NVLaserD 

0 

0 

0.3 

0.3 

0.6 

0.6 

1.2 

1.2 

a  (seeker  to  target  normal  angle  -  deg) 

0 

45 

0 

45 

0 

45 

0 

45 

Seeker  Signal  (x  10"5  W) 

sSalSeeker 

4.59 

3.25 

4.34 

3.07 

2.98 

2.11 

1.35 

0.96 

NVLaserD 

4.62 

3.26 

4.36 

3.08 

3.03 

2.14 

1.39 

0.98 

Difference  (%) 

-0.7 

-0.3 

-0.5 

-0.3 

-1.7 

-1.4 

-2.9 

-2.1 

Figure  17.  Comparison  of  2-a  and  1-ct  divergence. 

This  theory  is  supported  by  the  NVLaserD  beam  dimensions  at  the  target,  whose  standard 
deviations  were  twice  those  predicted  by  sSalSeeker.  Therefore,  to  compare  calculations  with 
the  NVLaserD  model,  the  sSalSeeker  5  values  were  made  double  those  of  NVLaserD.  After 
making  this  assumption,  the  models  agreed  within  3%  for  all  comparisons. 
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6.2  Geometry 
6.2.1  Projectile  Fly-In 

Initial  Conditions 


•  Target  Location  (tgeo)  =  (0,0,1.15) 

•  Designator  Location  (dgeo)  =  (-2000,0,1.15) 

•  Seeker  Location  (sgeo)  =  (-25000,0,25000) 

•  Seeker  Azimuth  (s_orientation.X)  =  0° 

•  Seeker  Elevation  (s_orientation.Y)  =  -45° 

•  Target  Rotation  (t  rotation)  =  0° 

Variables 

•  Seeker  Location  (sgeo)  =  (x,0,z)  -  varied  so  that  seeker  is  on  a  45°  descent  towards  target 

•  Visibility  (sal_weather. visibility)  =  clear,  hazy 


Figure  18  shows  the  results  of  a  virtual  fly-in  that  tested  how  much  the  laser  signal  increased  as 
the  range  to  target  decreased.  The  projectile  began  approximately  35  km  away  from  the  target 
(25  km  horizontally  and  25  km  vertically),  and  closed  in  on  the  target  at  a  45°  descent. 
Designator,  projectile,  and  target  were  all  aligned  along  the  x-axis  (no  side-to-side  movement). 
The  green  line  in  figure  1 8  represents  the  S/N  threshold  above  which  the  seeker  detects.  For  the 
fly-in  test,  the  seeker  detected  at  a  21.5-km  horizontal  range  for  clear  conditions  and  at  14.25-km 
horizontal  range  for  hazy  conditions. 


27 


Projectile  Fly-In 


Signal  (Clear)  Signal  (Hazy)  S/N  Threshold 


-25000  -23000  -21000  -19000  -17000  -15000  -13000  -11000  -9000  -7000  -5000 

Horizontal  Distance  from  Target  (m) 


Figure  18.  Projectile  fly-in  test. 

6.2.2  Off-Angle  Test 

Initial  Conditions 

•  Target  Location  (tgeo)  =  (0,0,1.15) 

•  Designator  Location  (dgeo)  =  (-2000,0,1.15) 

•  Seeker  Location  (sgeo)  =  (-2000,0,1.15) 

•  Seeker  Azimuth  (sorientation.X)  =  0° 

•  Seeker  Elevation  (s  orientation.Y)  =  0° 

•  Target  Rotation  (t  rotation)  =  0° 

Variables 

•  Seeker  Location  (sgeo)  =  varied  so  that  x  and  y  are  on  circle  with  radius  of  2000  m 

•  Seeker  Azimuth  (s  orientation.X)  =  varied  so  that  seeker  always  points  at  the  geometric 

center  of  the  target 
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The  seeker  was  positioned  at  a  series  of  points  on  a  circle  around  the  parallelepiped  target, 
always  pointing  at  the  target’s  geometric  center.  The  targets’  location  was  fixed  and  its  rotation 
was  0.  The  designator  remained  fixed,  and  aimed  along  the  +x  axis  to  designate  the  left  side  of 
the  target.  As  the  seeker's  y  deviated  from  0,  the  seeker's  off-angle  from  the  target  surface 
normal  (a,  see  figure  1)  increased  from  0°  to  90°  in  3-degree  increments,  and  the  received  power 
decreased  as  the  cosine  of  that  angle  (figure  19).  At  90°,  the  seeker  faced  parallel  to  the 
designated  surface  of  the  target,  at  the  boundary  of  the  hemisphere  containing  the  reflected  laser 
power.  Beyond  90°,  the  seeker  moved  outside  of  this  hemisphere,  and  the  received  power 
dropped  to  0.  Points  along  the  graph’s  horizontal  axis  are  seeker  positions  where  a  was  beyond 
90°. 


View  from  above  (z  points  out  of  page) 

♦  x 


Received  Power  vs  Seeker  Location 


Seeker  Location  (Y) 


Figure  19.  Off-angle  test. 
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6.2.3  Target  Rotation  Test 

Initial  Conditions 


•  Target  Location  (tgeo)  =  (0,0,1.15) 

•  Designator  Location  (dgeo)  =  (-2000,0,1.15) 

•  Seeker  Location  (sgeo)  =  (-5000,0,1.15) 

•  Seeker  Azimuth  (s  orientation.X)  =  0° 

•  Seeker  Elevation  (s  orientation.Y)  =  0° 

•  Target  Rotation  (trotation)  =  0° 

Variables 

•  Target  Rotation  (t  rotation)  =  0  to  360° 

•  Divergence  (d_divergence)  =  0.3,  0.6  mrad 

While  the  designator  and  seeker  stayed  at  a  fixed  distance  from  the  target,  the  target  rotation  was 
varied  from  0°  to  360°.  Figure  20  shows  how  the  power  varied  for  different  target  rotation 
angles,  with  the  rectangular  target  (6.4  x  2.3  m)  displaying  wider  variation  than  the  square  target 
(2.3  x  2.3  m).  For  6  =  0.3  mrad,  the  rectangular  target  underperformed  the  square  target  at  target 
rotation  angles  such  as  60°.  This  occurred  when  the  rectangle’s  long  surface  was  the  main  side 
receiving  power,  but  at  a  very  oblique  off-angle  (a)  from  the  seeker.  This  observation  poses  an 
interesting  question  for  the  designator  aim  point,  whether  it  is  best  to  point  at  the  target  center,  or 
at  the  side  most  nonnal  to  the  beam.  As  beam  divergence  increased,  the  larger  rectangular  target 
improved  its  perfonnance  relative  to  the  square  target,  as  the  beam  started  to  spill  off  of  the  sides 
of  the  smaller  square  target  (note  the  graph  for  5  =  0.6  mrad,  where  the  rectangular  target 
outperfonned  the  square  target  at  angles  of  0°  and  180°). 
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Rectangular  Ta  rget  Square  Target 
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Figure  20.  Target  rotation  test. 
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7.  Path  Forward 


In  summary,  the  model  calculates  the  power  distribution  across  a  SAL  seeker’s  IR  detector  given 
relevant  geometry,  laser  characteristics,  and  atmospheric  infonnation.  By  analyzing  this  power 
distribution,  the  model  returns  flight  guidance  information.  The  model’s  C++  class 
implementation,  sSalSeeker,  can  be  run  stand-alone,  and  is  also  easily  embeddable  into  larger 
smart-weapon  models. 

Recently,  researchers  have  investigated  the  ability  of  reduced-state  guidance  algorithms  to 
successfully  guide  a  munition  to  the  target  (10).  Upon  integration  into  a  larger  smart-weapon 
model,  this  model  could  provide  insight  into  another  simplified  guidance  scenario:  SAL-only 
guidance,  without  input  from  a  global-positioning  system  (GPS)  or  inertial  measurement  unit 
(IMU).  This  requires  the  ability  to  shoot  ballistically  into  a  -guidance  basket,”  after  which  the 
SAL  seeker  takes  over  (figure  21).  The  dimensions  of  the  basket  depend  largely  on  the  seeker’s 
FOV,  the  rate  at  which  the  laser  beam  attenuates,  and  the  maneuver  authority  of  the  projectile. 


Figure  21.  Extending  the  guidance  basket. 
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The  following  are  suggestions  for  further  development  of  the  model: 


•  Target  geometry 

Current  -  Rectangular  parallelepiped 

Future  -  Detailed  target  geometries  including  windows  and  tires  would  permit  study  of 
whether  there  is  an  ideal  location  to  designate  for  different  targets. 

Difficulty  -  Low 

•  Terrain 

Current  -  Flat  terrain,  which  makes  for  an  unrealistically-favorable  situation. 

Future  -  Variable-height  terrain,  with  vegetation  and  man-made  features 
Difficulty  -  Medium 

•  Target  reflection 

Current  -  Assumes  perfect  Lambertian  surface  reflection 

Future  -  Combination  of  Lambertian  and  specular  reflection,  depending  on  the  surface. 
Reflection  would  then  depend  on  the  angle  of  incidence  from  the  designator.  For 
example,  as  target  geometries  become  more  developed,  windows  may  incorporate  more 
specular  reflection  than  matte-finished  doors. 

Difficulty  -  Medium 

•  Background  noise 

Current  -  Input  parameter 

Future  -  Calculation  of  the  background  noise  from  the  sun  and  earth  entering  the 
seeker 

Difficulty  -  Medium 
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Appendix  A.  Sample  Input  and  Output 


This  appendix  is  in  its  original  form,  without  editorial  change. 
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Figure  A-l.  Sample  input  file. 
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-ia|x| 

SAL  SEEKER  CLASS  IMPLEMENTATION 

- 

developed  for  AWCB  SWEEPM  (LSS;  2011) 

.  INPUT  PARAMETERS  . 

WEATHER 
Latitude, 

Season, 

Visibility, 

SEEKER 
Seeker  FOV  (Half  Angle), 
Seeker  Aperture  Diameter, 


***  DESIGNATOR 
Laser  Type, 

Laser  Divergence, 
Ray  Count, 
Designator  Error 
Designator  Error 


83: 


*  TARGET 

Size  (Length;  width;  Height), 
Reflectivity, 

Left  Side, 

Front, 

Right  Side, 

Sack, 

Top, 

Ground , 


Designator  to  Target  Distance, 
Seeker  to  Target  Distance, 

Seeker  Location  (XYZ  -  m) , 
Designator  Location, 

Target  Location, 

Seeker  Heading, 

Designator  Heading, 

Target  Rotation, 

Seeker  Azimuth, 

Seeker  Elevation, 


Surface  Most  Hit, 
Left  Side, 

Front , 

Right  Side, 

Sack, 

Top, 

Ground, 

Di rect IntoSeeker , 
<iss. 


SURFACE  HITS 


Ql  (+  Pitch;  +  Yaw) 
Q2  (+  Pitch;  -  Yaw) 
Q3  C-  Pitch;  +  Yaw) 
Q4  (-  Pitch;  -  Yaw) 


SEEKER  QUADRANT  HITS 


Designator  Laser  Signal, 
Seeker  Laser  Signal, 
Background  Noise, 

1/N  Threshold, 

:/N  Ql, 

S/N  Q2, 

S/N  Q3, 

S/N  Q4 , 

S/N  Total , 


WWhK  LOSS 


0  =  No;  1  =  Yes 
Target  Encountered, 
Target  Detected, 


ENCOUNTER/DETECT 


QUADRANT  SIGNALS 
+  Pitch 


Mid-lat 

Summer 

Clear 


4.50  degrees 
60.00  mm 


Nd:Glass 
0. 30  mrad 
10000 
0.10  mrad 
0.10  mrad 


6.4,  2.3,  2.3m 

0.40 

0.40 

0.40 

0.40 

0.40 

0.40 


2.00  km 
5.00  km 

-5000,  0,  1 
-2000,  0,  1 
0.  0,  1 

<1.00,  0.00,  0.00> 
<1.00,  -0.00,  -0.00> 
0  degrees 
0  degrees 
0  degrees 


1284 

2101 

2524 

4091 


5. 33e+006  W 

4. 59e-005  W 

5.48e-007  W 

7.00 

43.07 

70.48 

84.63 

136.93 

83.78 


‘J 


Ql 

12.85  % 

Q2 

21.03  * 

Q3 

25.26  % 

Q4 

40.86  % 

-  Pitch 

MANEUVER  GUIDANCE 


Pitch  Command, 

Yaw  Command, 

Signal  Center:  Theta, 
Signal  Center:  Magnitude, 


-0.32 

-0.24 

306.42 

0.28 


Press  any  key  to  continue 


A 

jJ 


Figure  A-2.  Sample  output  file. 
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Intentionally  left  blank. 
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This  appendix  is  in  its  original  form,  without  editorial  change. 
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s_sal_seeker_class.Ii 

/* 

SAL  Seeker  Class  . h 

Luke  Strohm,  Army  Research  Laboratory 
3-14-11 
*/ 


/******************************************************************************/ 
#ifndef  S_SAL_SEEKER_CLASS_H_ 

#def ine  S_SAL_SEEKER_CLASS_H_ 

!  • ************************************************************************ 
#include  <iostream> 

#include  <iomanip> 

#include  <io.h> 

#include  <fstream> 

#include  <string> 

#include  <cstdlib> 

#include  <windows.h> 

#include  <cmath> 

j  ' ************************************************************************ 

#include  "utilities  .  h"  // . MKA  Utility 

#include  "statistics  .h"  // . MKA  Statistics 

#include  "y_format_namespace  .h"  // . RJY  Formatting 

using  namespace  std; 
using  namespace  yFormat; 

/ 

II 

st 


★  -*•-*•**-*•  j 

Library 

Library 

Library 


Define  Structs 
ruct  weather { 

int  season;  // . 

// . 

int  latitude;  //... 

/  / ... 
n... 

int  visibility;  //. 

II. 


. . 0  =  summer  (3/22-9/21) 

. 1  =  winter  (9/22-3/21) 

. . 0  =  tropics  (<23.5  deg  +-) 

. 1  =  mid-latitudes  (23.5-50  deg  +-) 

. . 2  =  sub-arctic  (50-70  deg  +-) 

. . 0  =  clear  (23  km  visibility) 

. 1  =  hazy  (5  km  visibility) 


}  ; 

st 


ruct  rayhit{ 

Point3D  hit;  // . X,  Y,  Z  Coords  of  where  ray  intersects  surface 

int  surf;  // . Surface  point  hits 

double  power;  // . Power  of  ray 

bool  obscured;  //.Check  to  see  if  ground  hit  is  obscured  from  seeker  by  target 


}  ; 

st 


ruct  quadsignal{ 
double  ql;  II... 
double  q2;  // . . . 
double  q3;  //.... 
double  q4 ;  // . . . 


//. 


.Contains  info  on  the  four  quads  of  quad-band  detector 

. Positive  Pitch,  Positive  Yaw 

. Positive  Pitch,  Negative  Yaw 

. Negative  Pitch,  Positive  Yaw 

. Negative  Pitch,  Negative  Yaw 


}  ; 

st 


ruct  guidancesignal {  // . Actuator  guidance  (-1/  +  1  is  maximum  maneuver) 

double  pitch;  // . Pitching  (Lifting!)  Multiplier  (-1  to  1) 

double  yaw;  // . Yawing  (Left  Turn!)  Multiplier  (-1  to  1) 

double  theta;  // . Polar  Direction  of  Signal  Center 

// . 0-2pi  rad,  0  =  -yaw  dir,  rotate  CCW 

double  mag;  // . Polar  Magnitude  of  Direction  Vector  (0-1) 

ruct  saloutputs{ 
bool  encounter; 
bool  detect; 
quadsignal  sal_signals; 
guidancesignal  actuator  signals; 

} ; 

j^^^^^^^^^^-k-k-k-k'k-k-k'k-k-k-k-k-k-k-k'k-k-k-k-k'k-k-k'k-k-k-k-k'k-k-k'k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k'k-k-k-k-k-k-k-k'k-k-k'k-k-k-k-kl 

class  sSalSeeker{ 


}  ; 

st 
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public : //************************ - PUBLIC - *********************************** 


//****************************0UTSIDE  CLASSES* ************************  ******** 

Utilities  UTILITY;  // . Instantiate  Utility  Class  (MKA) 

Statistics  STAT;  // . Instantiate  Statistics  Class  (MKA) 

/ /********************************OUTPUT** ******************************  ****** 

/ / 0_FUNCTI0NS - 

sSalSeeker ( ) ; 

saloutputs  S_SALSEEKER (Point3D  d_geo,  Point3D  t_geo,  //..Runs  all  calculations 
Point3D  s_geo,  Point2D  s_orientation, // . and  outputs  to  struct 
double  t_rotation,  int  seed) ; 
saloutputs  programoutput ; 

void  PrintOutput ( string  filename  =  string  format  =  ".txt");  //  SummaryFile 

/ / 0_DI SPLAY - 

bool  printtoscreen; 
bool  printtofile; 

int  surface[8];  / / . 0  =  Left  Side,  1  =  Front,  2  =  Right  Side,  3  =  Back,  4  =  Top 

// . 5  =  Ground  (underspill),  6  =  DirectlntoSeeker ,  7  =  Miss 

string  s_desc[8];  //.Create  array  of  target  surface  names  (ref  above  surfaces) 
int  mainhitsurface ;  // . Surface  that  got  hit  the  most 

//**** ******  ************  ******** parameters** **************************  ******** 


/ 


//PJ3EEKER - 

double  s_fov;  // . deg  (half  angle  of  cone  swept  out) 

double  s_aperture_diameter ;  // . m 

double  s_background;  // . background  noise  multiplier  W/mA2/sr 

double  s_signalnoise_threshold;  //ratio  of  signal  to  noise  necessary  to  detect 

double  s_eff iciency;  // . efficiency  factor  accounting  for  various  losses 

/ / P_DE  S I GNAT  OR - 

int  d_lasertype;  // . 0  =  1.06  urn  (Nd:  Glass) 

// . 1  =  1.536  urn  (Er:Glass) 

int  d_raycount;  // . Number  of  rays  to  divide  laser  pulse  into 

double  d_divergence;  // . rad 

double  d_pulse_energy;  // . J 

double  d_pulse_f requency ;  // . Hz 

double  d_pulse_duration;  // . s 

double  d_eff iciency;  // . efficiency  factor  accounting  for  various  losses 

double  d_h_error,  d_v_error;  //.... designator  error  (hor  and  vert,  rad,  stdev) 

/  /  P_TARGET  - - — - 

Point3D  t_size;  // . length,  width,  height:  m 

double  t_reflect [6] ;  //  target  reflectivity  for  each  surface  (ref  surface  des) 

//PJ3THER - 

int  seed;  // . seed  number  for  random  number  generator 

weather  sal_weather;  // . struct  containing  weather  information 

double  s_noise;  // . W 

★  **★★**★★**★★**★★**★★***★**★★**★★**★★**★★**★★**★★**★★***★**★★**★★**★★**★★**★•*■*  j 


private: //************************ - PRIVATE - ********************************* 

//**** ****************************VARIABLES* ************************  ********** 

/* 

Rotations  use  RIGHT-HANDED  coordinate  system 
+  z  +y 

I  / 

!  / 

I  / 

I  / 

I  / 

-------  +x 

XYZ  Zero  Points: 
x  =  0  -  user-defined 
y  =  0  -  user-defined 
z  =  0  -  earth' s  surface 
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Azimuth  is  angle  between  +x  direction  and 
projection  of  flight  body  onto  x-y  plane 
+  Azimuth  =  CCW  rotation  as  seen  from  above 


Elevation  is  angle  between  flight  body  and  x-y  plane 
+  Elevation  =  Flight  body  above  x-y  plane 

*/ 

Point3D  d_geo;  // . Designator  location  projected  onto  Earth-faced 

// .  x/y/z  coords  (m)  (See  further  description  above) 

Point3D  t_geo;  // . Target  location  projected  "  "  (m) 

Point3D  s_geo;  // . Target  Location  projected  "  "  (m) 

Point2D  s_orientation;  // . Azimuth/Elevation  of  projectile  (deg) 

// . (See  further  description  above) 

double  t_rotation;  // . deg  (0  degrees  is  the  target  maneuvering 

// . in  the  +y  direction,  +  deg  rotates 

// . target  clockwise  as  seen  from  above) 


-k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 


// 

//F  MAIN- 


FUNCTIONS 


•k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k-k 


void  SetDefaults  ( ) ; 
void  LaserSpot  //... 
(Point3D  pi,  Point3D 

// . 

// . 

// . 

// . 

// . 

void  SignalLoss ( ) ; 


//. 


. Sets  default  values  of  initial  parameters 

...Populates  p  vector  with  power,  hit  points,  and  surfaces 
p2,  Point3D  p3,  Point2D  s  or,  double  tr) ; 

. pi  =  Designator  Center  (x,y, z) 

. p2  =  Target  Center  (x,y, z) 

. p3  =  Seeker  Center  (x,y,z) 

. s_or  =  seeker  orientation  (azimuth,  elevation) 

. tr  =  Target  Rotation 

// . Calculates  signal  loss  through  attenuation, 

// . reflection,  and  designator/receiver  losses 

//FJJTILITY - 

//  Calculates  attenuation  through  atmosphere  via  lookup  table 

double  Attenuation (double  currentsig, double  traveldistance, double  startheight, 
double  endheight,  weather  sal_weather,  int  lasertype) ; 
Point3D  SpherCartConv (Point2D  pi);//  Converts  azimuth/elevation  to  unit  vector 
//  Creates  unit  vector  pointing  from  pi  to  p2  (normalized  by  default) 

Point3D  VectorPointing (Point3D  pi,  Point3D  p2,  bool  norm  =  true) ; 

//  Rotates  vector  by  angle  a  (yaw)  and  b  (pitch) 

Point3D  VectorRotation (Point3D  pi,  double  a,  double  b,  int  vect  =  0) ; 

double  VectorAngle (Point3D  pi,  Point3D  p2) ;  // . Angle  between  two  vectors 

double  DotProduct (Point3D  pi,  Point3D  p2);  // . Dot  product  of  two  vectors 

void  InverseMatrix (double  m[3] [3], double  mi [3] [3], double  &det) ; //Inversematrix 
//******************************  CALCULAT IONS********************************** 

//C_TARGET_GEOMETRY - 

vector  <Point3D>  s,  n;  // . s  =  Target  surface  coordinates 

// . n  =  Surface  normal  vectors 

vector  <rayhit>  r;  // . Struct  containing  ray  hit  information 

/ / C_DIRECTION_VECTORS - 

Point3D  d_heading_vector ;  // . Designator  pointing  vector 

Point3D  s_heading_vector ;  // . Seeker  pointing  vector 

/  /  C__ENERGY_LOS  S - 

double  s_aperture_area;  // . mA2 

double  s  t  multiplier;  //..Ratio  of  power  reaching  seeker  aperture  from  target 

/ / C_SEEKER_S I GNAL - 

double  s_power;  // . W 

int  s_hits[4];  // . number  of  ray  hits  received  in  each  quadrant 

quadsignal  s_signal;  // . Total  signal  received  in  each  quadrant 

quadsignal  s_signalnoise;  // . Signal/Noise  Ratio  for  each  quadrant 

guidancesignal  smaneuver;  //..Pitching  and  yawing  maneuver  guidance  (-1  to  1) 

//CJ3ETECT - - - 

bool  s  encounter;  // . True  if  target  is  within  seeker's  FOV 

bool  s_detect;  // . True  if  laser  S/N  >  threshold 

}./****************************************************************************/ 

#endif 
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/  ' k-k-k-k'k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k'k-k-k'k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k'k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k  , 


s_sal_seeker_class.cpp 

/* 

SAL  Seeker  Class  . cpp 

Luke  Strohm,  Army  Research  Laboratory 
3-14-11 
*/ 


/******************************************************************************/ 
# include  " s_sal_seeker_class . h" 

I****************************************************************************** / 

sSal Seeker : : sSal Seeker  ( )  { 

SetDefaults ( ) ; 

}/*****************************************************************************/ 
void  sSalSeeker :: SetDefaults  ()  { 

//  Output  Parameters 
printtoscreen  =  false; 
printtofile  =  false; 

//  Weather  Parameters 

sal_weather . latitude  =  1;  // . 0  =  tropics  (<23.5  deg  +-) 

// . 1  =  mid-latitudes  (23.5-50  deg  +-) 

// . 2  =  sub-arctic  (50-70  deg  +-) 

sal_weather . season  =  0;  // . 0  =  summer  (3/22-9/21) 

// . 1  =  winter  (9/22-3/21) 

sal_weather . visibility  =  0;  // . 0  =  clear  (23  km),  1  =  hazy  (5  km) 

//  Seeker  Parameters 

s_fov  =  4.5;  s_fov  *=  UTILITY . DEG2RAD;  //...deg  (half  angle  of  cone  swept  out) 

s_aperture_diameter  =  60E-3;  // . m  (based  on  81mm  mortar) 

s_background  =  0.010;  // . W/mA2/sr 

s_signalnoise_threshold  =  7.0; 

s_efficiency  =  0.95;  // . efficiency  factor  encompassing  various  losses 

//  Designator  (Scout)  Parameters 

d_lasertype  =  0;  // . 0  =  1.06  um  (Nd;  Glass) 

// . 1  =  1.536  um  (Er:Glass) 

d_raycount  =  10000;  // . number  of  rays  in  laser  beam 

d_divergence  =  3E-4;  // . rad 

d_pulse_energy  =  80E-3;  // . J 

d_pulse_f requency  =  10.0;  // . Hz 

d_pulse_duration  =  15E-9;  // . s 

d_efficiency  =  0.95;  // . efficiency  factor  encompassing  various  losses 

d_h_error  =  IE-4;  // . horizontal  designator  error  (rad,  stdev) 

d  v  error  =  IE-4;  // . vertical  designator  error  (rad,  stdev) 

//  Target  Parameters 


t 

si 

ze .  X 

=  6 

4, 

//. 

t 

si 

ze .  Y 

=  2 

3, 

//. 

t 

si 

ze .  Z 

=  2 

3, 

//. 

t 

re 

fleet 

[0] 

= 

0 

.4; 

// 

t 

re 

fleet 

[1] 

= 

0 

.4; 

// 

t 

re 

fleet 

[2] 

= 

0 

.4; 

// 

t 

re 

fleet 

[3] 

= 

0 

.4; 

// 

t 

re 

fleet 

[4] 

= 

0 

.4; 

// 

t 

re 

fleet 

[5] 

= 

0 

.4; 

// 

. m 

. m 

. m 

.  .  left  side 

. front 

.  right  side 

. back 

. top 

. ground 


j/*****************************************************************************/ 

saloutputs  sSalSeeker : : S_SALSEEKER (Point3D  d_geo,  Point3D  t_geo,  Point3D  s_geo, 

Point2D  s_orientation,  double  t_rotation,  int  seed) { 
//  Set  class  variables  equal  to  function  inputs 
(*this) ,d_geo  =  d_geo; 

(*this) . t_geo  =  t_geo; 

(*this) . s_geo  =  s_geo; 

( *this )  . s_orientation  =  s_orientation; 

(*this).t  rotation  =  t  rotation; 


//  Initialize  random  number  generator 
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STAT .MIRandGenlnit (seed)  ; 


//  Call  Calculation  Functions 

LaserSpot (d_geo,  t_geo,  s_geo,  s_orientation,  t_rotation) ; 

SignalLoss  ( ) ; 

programoutput . encounter  =  s_encounter; 
programoutput . detect  =  s_detect; 
programoutput . sal_signals  =  s  signal ; 
programoutput . actuator  signals  =  s_maneuver; 

return  programoutput; 

}/*****************************************************************************/ 
void  sSalSeeker : : LaserSpot 

(Point3D  pi,  Point3D  p2,  Point3D  p3,  Point2D  s_or,  double  tr) { 

//  Define  local  variables 

vector  <Point3D>  u(7),  v(7);  //..u  and  v  (vectors  from  surface  center  to  edge) 

// . SIDES:  u  (length  dir),  v  (height  dir) 

// . FRONT,  BACK:  u  (width  dir),  v  (height  dir) 

// . TOP,  BOTTOM:  u  (length  dir),  v  (width  dir) 

Point3D  r_ind;  // . individual  ray  vector  from  designator  to  target 

double  m[3] [3] ,  mi [3] [3],  dl[3],  tuv[3]; 

double  det  =  0.0;  // . determinant 

double  r_mult  =  0.0;  // . multiplier  to  get  from  designator  to  surface  hit 

int  hitcount  =  0;  // . the  number  of  hits  on  the  surface  most  hit 

s_desc[0]  =  "Left  Side";  s_desc[l]  =  "Front";  s_desc[2]  =  "Right  Side"; 
s_desc[3]  =  "Back";  s_desc[4]  =  "Top";  s_desc[5]  =  "Ground"; 
s_desc[6]  =  "DirectToSeeker" ;  s_desc[7]  =  "Miss"; 

//  Resize  vectors  s,  n,  and  r 

s.resize(7);  n.resize(7);  r . resize (d_raycount) ; 

//  Initialize  surface  hit  counts  to  0 

surface[0]  =  surfacefl]  =  surface[2]  =  surface[3]  =  surface[4] 

=  surface[5]  =  surface[6]  =  surface[7]  =  0; 

//  Target  Rotation  -  only  valid  0  to  2PI  rad 
if  ( (tr  <  0)  ||  (tr  >  2*M_PI ) ) { 

cerr  <<  "Invalid  target  angle"  «  endl; 
system ("PAUSE") ; 
exit  (1 ) ; 

} 


//  DEFINE  TARGET  SURFACE  COORDINATES  AND  VECTORS***************************** 

//  Calculate  Heading  (Pointing)  Vector  for  Seeker 
s_heading_vector  =  SpherCartConv ( s_or ) ; 


// 

Target  Surface  Coordinates 

// 

Side  0  (Left  Side) 

s[0] 

.X  -  p2.X  -  (t_size . Y/2 ) 

* 

cos  (tr) 

s[0] 

.Y  =  p2 . Y  +  (t  size. Y/2) 

sin  (tr) 

s  [0] 

. Z  =  p2 . Z; 

// 

Side  1  (Front) 

S[l] 

.X  =  p2 .X  +  (t  size . X/2 ) 

* 

sin  (tr) 

S[l] 

.Y  =  p2 . Y  +  (t  size. X/2) 

: k 

cos  (tr) 

S[l] 

| . Z  =  p2 . Z; 

// 

Side  2  (Right  Side) 

s[2] 

.X  =  p2.X  +  (t  size. Y/2) 

* 

cos  (tr) 

s  [  2  ] 

I . Y  —  p2 . Y  -  (t  size. Y/2) 

Vc 

sin  (tr ) 

s[2] 

. Z  =  p2 . Z; 

44 


sin  (tr) ; 
cos  (tr)  ; 


//  Side 
s [3] .X  = 
s  [ 3 ]  .Y  = 
s  [  3  ]  .  Z  = 
//  Side 
s [4] .X  = 
s [4] .Y  = 
s [4] .Z  = 
//  Side 
s [5] .X  = 
s [5] .Y  = 
s  [5]  .  Z  = 
//  Side 
s [6] .X  = 
s [6] .Y  = 
s  [  6  ]  .  Z  = 


3  (Back) 

p2.X  -  (t_size.X/2) 
p2.Y  -  (t_size.X/2) 
p2  .  Z; 

4  (Top) 
p2  .X; 
p2  .  Y; 

p2 . Z  +  (t_size.Z/2) 

5  (Ground) 
p2  .X; 

p2  .  Y; 


0; 

6  ( DirectToSeeker ) 
p3 .  X; 
p3.  Y; 
p3 .  Z ; 


* 

* 


// 

// 

n  [  0] 
n  [  0] 
n  [  0  ] 
// 

n  [  1] 
n  [  1  ] 
n  [  1  ] 

// 

n  [  2  ] 
n  [2] 
n  [2] 
// 

n  [  3] 
n  [3] 
n  [3] 
// 

n  [4] 
n  [4] 
n  [4] 
// 

n  [5] 
n  [5] 
n  [5] 
// 

n  [  6] 
n  [  6  ] 
n  [  6] 


Normal  Vectors  Pointing  out  of  each  surface 
Side  0  (Left  Side) 

.X  =  -cos (tr) ; 

. Y  =  sin  (tr )  ; 

.  Z  =  0; 

Side  1  (Front) 

.X  =  sin  (tr)  ; 

. Y  =  cos  (tr ) ; 

.Z  =  0; 

Side  2  (Right  Side) 

.X  =  cos  (tr) ; 

. Y  =  -sin (tr) ; 

.Z  =  0; 

Side  3  (Back) 

.X  =  -sin  (tr) ; 

-cos  (tr ) ; 

0; 


.Y  = 
.Z  = 
Side 
.X  = 
.Y  = 
.Z  = 
Side 
.X  = 
.Y  = 
.Z  = 
Side 
.X  = 
.Y  = 
.Z  = 


4  (Top) 

0; 

0; 

1; 

5  (Ground) 

0; 

0; 

1; 

6  (DirectToSeeker) 
s_heading_vector . X; 
s_heading_vector . Y; 
s  heading_vector . Z ; 


// 

Surface-defining  vectors  (u 

// 

Side 

0  (Left  Side) 

u  [  0  ] 

.X  = 

n[l].X  *  t  size.X/2; 

u  [  0  ] 

.Y  = 

n[l].Y  *  t  size.X/2; 

u  [  0  ] 

.Z  = 

0; 

v[0] 

.X  = 

0; 

v  [  0  ] 

.Y  = 

0; 

v  [  0  ] 

.Z  = 

t  size.Z/2; 

// 

Side 

1  (Front) 

U[l] 

.X  = 

n[2].X  *  t  size.Y/2; 

U[l] 

.Y  = 

n[2].Y  *  t  size.Y/2; 

U[l] 

.Z  = 

0; 

v[l] 

.X  = 

0; 

v[l] 

.Y  = 

0; 

v[l] 

.Z  = 

t  size.Z/2; 

// 

Side 

2  (Right  Side) 

and  v) 
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u  [  2  ]  .X 
u  [  2  ]  .Y 
u  [  2  ]  .Z 
v  [  2  ]  .X 
v  [  2  ]  .Y 
v  [  2  ]  .Z 
//  Side  3 
u[3] .X  =  n 
u  [  3  ]  .Y  = 
u  [  3  ]  .Z  = 
v  [  3  ]  .X  = 
v  [  3  ]  .Y  = 
v  [  3  ]  .Z  = 

//  Side 
u  [  4  ]  .X  = 
u  [  4  ]  .Y  = 
u  [  4  ]  .  Z  = 
v  [  4  ]  .X  = 
v  [  4  ]  .Y  = 
v  [  4  ]  .  Z  = 

//  Side 
u [5] .X  = 
u [5] .Y  = 
u  [5]  . Z  = 
v  [  5]  .  X  = 
v  [  5  ]  .Y  = 
v  [  5]  . Z  = 

//  Side 
//  Define 
Point3D  vl 
Point3D  v2 
Point3D  v3 


[3].X  *  t_size.X/2; 
[3 ] . Y  *  t  size . X/2 ; 


_size .  Z/2 ; 

(Back) 

[0].X  *  t_size.Y/2; 
[ 0 ] . Y  *  t  size.Y/2; 


_size .  Z/2 ; 

(Top) 

[1] .X  *  t_size .  X/2 ; 

[1]  .Y  *  t  size. X/2; 

[2]  .X  *  t  size.Y/2; 

[2].Y  *  t  size.Y/2; 

(Ground) 

[1] .X  *  t_size .  X/2 ; 

[1]  .Y  *  t  size. X/2; 

[2]  .X  *  t  size.Y/2; 

[2] .Y  *  t  size.Y/2; 

( Direct To Seeker ) 

vectors  that  characterize  seeker  pointing 
=  s_heading_vector ; 

=  VectorRotation (vl, 0, 0, 1) ; 

=  VectorRotation (vl, 0,0,2); 


u [ 6 ] . X  =  v2 . X 
u [ 6 ] . Y  =  v2 . Y 
u [ 6 ] . Z  =  v2 . Z 
v  [  6 ]  . X  =  v3 . X 
v  [  6 ]  . Y  =  v3 . Y 
v  [  6 ]  . Z  =  v3 . Z 


*  (s_aperture_diameter  /  2); 

*  (s_aperture_diameter  /  2); 

*  (s_aperture_diameter  /  2); 

*  (s_aperture_diameter  /  2); 

*  (s_aperture_diameter  /  2); 

*  (s_aperture_diameter  /  2); 


//  DIVIDE  LASER  BEAM  INTO  INDIVIDUAL  RAYS************************************ 


//  Calculate  Pointing  Vector  for  Designator  (with  Error) 
d_heading_vector  =  VectorRotation (VectorPointing (d_geo,  t_geo) , 

STAT .Normal (0, d_h_error)  , 
STAT .Normal (0, d  v  error) )  ; 

for  (int  i=0 ; i<d_raycount;  i  +  +  )  { 

/*  Normally  distribute  rays  according  to  divergence 

Divergence  for  Gaussian  beams  assumed  to  be  2  st  deviations 
Beam  edge  when  intensity  drops  to  I  =  Io/eA2  */ 

r_ind  =  VectorRotation (d_heading_vector, 

STAT. Normal (0, d_divergence/2) , 

STAT. Normal (0, d_divergence/2) ) ; 


//  Set  default  values 

r_mult  =  1E10;  // . Any  hits  should  produce  r_mults  under  1E10  m 

int  tempsurface  =  7;  // . Surface  7  means  that  ray  did  not  hit  anything 

//  FIND  WHAT  SURFACE  RAY  HITS********************************************** 
for  (int  j=0; j<7; j++) { 

//  Calculate  m 
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m [ 0 ] [0]  =  -r_ind.X;  m[0] [1]  =  u[j] .X;  m[0] [2]  =  v[j] .X; 

m [1] [0]  =  -r_ind.Y;  m[l] [1]  =  u[j] .Y;  m[l] [2]  =  v[j] .Y; 

m [ 2 ] [0]  =  -r_ind.Z;  m[2] [1]  =  u[j] .Z;  m[2] [2]  =  v[j] .Z; 

//  Calculate  mi  and  determinant 
Inver seMatrix (m, mi , det ) ; 

//  Calculate  dl 
dl [0]  =  pi .X  -  s [j ] .X; 
dl [1]  =  pl.Y  -  s [j] .Y; 
dl [2 ]  =  pl.Z  -  s  [  j ] .Z; 

//  Calculate  TUV  vector  (T  =  multiplier,  U  =  diml ,  V  =  dim2) 
tuv [ 0 ]  =  mi [0] [0] *dl [0]  +  mi [0 ] [ 1 ] *dl [ 1 ]  +  mi [0] [2] *dl [2 ] ; 

tuv [ 1 ]  =  mi [1]  [0] *dl [0]  +  mi [1]  [1] *dl  [1]  +  mi [1] [2] *dl [2] ; 

tuv [ 2 ]  =  mi  [2]  [0] *dl [0]  +  mi [2 ]  [1] *dl  [1]  +  mi [2] [2 ] *dl [2 ] ; 


} 


//  Determine  if  ray  intersects  plane 

if  ( (fabs (det) <lE-20)  ||  tuv[0]  <=  IE-20)  continue; 

//  Determine  if  ray  intersects  surface  on  plane 

if  (j  !=  5){  // . Ground  plane  is  infinite 

//  Check  intersection  straight  to  seeker 


if  (j  ==  6) { 

double  dist  =  sqrt (tuv [1 ] *tuv [1]  +  tuv [2 ] *tuv [2 ] ) ; 
double  fov_check  =  VectorAngle (s_heading_vector, 

VectorPointing ( s_geo,  d_geo) ) ; 

if  ((dist  >  s_aperture_diameter/2)  | |  (fov^check  >  s_fov) )  continue; 

} 

//  Rectangular  Target 

else  if  ((tuv[l]  <  -1)  |  |  (tuvfl]  >1)  I  I 

(tuv  [2]  <  -1)  ||  (tuv  [2]  >  1))  continue; 

} 

//  Look  for  lowest  r_mult 
if  ( tuv [ 0 ]  <  r_mult) { 
r_mult  =  tuv[0]; 
tempsurface  =  j ; 


} 


//  RECORD  HIT  POINT  AND  SURFACE  MOST  HIT*********************************** 

//  Add  ray  hit  point  to  array  of  hit  points  (p) 

if  (tempsurface  ==  6)  {  // . Ray  goes  directly  into  seeker 

r[i].hit.X  =  p3.X; 

r[i].hit.Y  =  p3.Y; 

r [ i]  . hit . Z  =  p 3 . Z ; 

} 

else  if  (tempsurface  ==  7){  // . Ray  hits  nothing 

r [ i ] . hit .X  =  pl.X; 
r [ i ] . hit . Y  =  pl.Y; 
r[i] .hit . Z  =  pl.Z; 

} 

else{  // . Ray  reflecting  off  of  surface 

r[i] .hit.X  =  pl.X  +  r_ind.X  *  r_mult; 

r[i] .hit.Y  =  pl.Y  +  r_ind.Y  *  r_mult; 

r[i] .hit.Z  =  pl.Z  +  r_ind.Z  *  r_mult; 


r[i]  .surf  =  tempsurface;  // . Record  surface  hit  by  ray 

surface  [tempsurface] ++;  // . Tally  hits  on  each  surface 


//  Check  ground  hits  to  see  if  they  are  obscured  from  seeker  by  target 
r [i] .obscured  =  false; 
if  (tempsurface  ==  5) 

{ 

Point3D  obs  =  VectorPointing (r [i ]. hit,  s_geo) ; 
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double  obsc_mult  =  1E10; 
int  obsc  surface  =  5; 


} 


for  (int  k=0 ; k<7 ; k++ ) { 

//  Calculate  m 

m [ 0 ] [0]  =  -obs.X;  m[0] [1]  =  u[k] .X;  m[0] [2] 

m [ 1 ] [0]  =  -obs.Y;  m[l] [1]  =  u[k] .Y;  m[l] [2] 

m [ 2 ] [0]  =  -obs.Z;  m[2] [1]  =  u[k] .Z;  m[2] [2] 

//  Calculate  mi  and  determinant 


v [ k] .X; 
v [ k] .Y; 
v  [  k]  .  Z; 


InverseMatrix (m, mi , det ) ; 

//  Calculate  dl 
dl  [0]  =  r  [i]  .hit . X  -  s  [k]  .X; 
dl [1 ]  =  r [i] .hit . Y  -  s[k].Y; 
dl [2 ]  =  r [ i ] .hit . Z  -  s[k] .Z; 

//  Calculate  TUV  vector  (T  =  multiplier,  U  =  diml ,  V  =  dim2 ) 
tuv [ 0 ]  =  mi [0] [0] *dl [0]  +  mi [0 ] [ 1 ] *dl [ 1 ]  +  mi [0] [2 ] *dl [2] ; 

tuv [ 1 ]  =  mi  [1]  [0] *dl [0]  +  mi [1 ]  [ 1 ] *dl  [  1 ]  +  mi  [1]  [2 ] *dl [2 ]  ; 

tuv [ 2 ]  =  mi [2]  [0] *dl [0]  +  mi [2 ]  [ 1 ] *dl  [  1 ]  +  mi [2 ]  [2 ] *dl [2 ]  ; 

//  Determine  if  ray  intersects  plane 

if  ( (fabs (det) <lE-20)  II  tuv[0]  <=  IE-20)  continue; 


//  Determine  if  ray  intersects  surface  on  plane 
if  ( ( tuv [ 1 ]  <  -1)  ||  ( tuv [ 1 ]  >1)|| 

(tuv[2]  <  -1)  ||  (tuv[2]  >  1))  continue; 


//  Look  for  lowest  obsc_mult 
if  (tuv[0]  <  obsc_mult) { 
obsc_mult  =  tuv[0]; 
obsc_surface  =  k;  } 

} 

//  Determine  if  ray  hit  is  obscured  from  seeker  by  target 
if  (obsc  surface  <  5)  r[i] .obscured  =  true; 


//  Determine  "mainhitsurf ace" 
for  (int  1=0;1<8;1++) 


if  (surface [1]  >  hitcount) 


hitcount  =  surface [1]; 
mainhitsurface  =  1; 

} 

} 


}/*****************************************************************************/ 
void  sSalSeeker ; : SignalLoss  ( ) { 

//  Initialize  Quadrants 

s_signal.ql  =  s_signal.q2  =  s_signal.q3  =  s_signal.q4  =  0; 
s_hits[0]  =  s_hits[l]  =  s_hits[2]  =  s_hits[3]  =  0; 

//  Initialize  S_Encounter  and  S_Detect  to  False  (Default) 
s_encounter  =  s_detect  =  false; 

//  Define  vectors 

Point3D  d_t_vector;  // . Distance  between  des  and  target  (km) 

Point3D  t_s_vector;  // . Vector  from  target  dot  to  seeker 

Point3D  s_t_vector;  // . Vector  from  seeker  to  target  dot 

Point3D  t_normal_vector ;  //....Vector  pointing  out  of  the  ilium  target  surface 

//  Define  Distances/Rotations/Angles 

double  d_t_distance;  //...Straight-line  distance  between  desig  and  target  (km) 
double  s_t_distance;  //..Straight-line  distance  between  seeker  and  target  (km) 

double  s_t_normal_angle;  // . Angle  betwn  seeker  pointing  and  target  normal 

double  s_t_correct_angle;  // . Angle  that  seeker  is  off  from  ray  hit 

double  r_correct_angle ;  // . Individual  FOV  check  for  each  ray 
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//  ENCOUNTER:  Check  to  see  if  target  is  encountered  by  seeker 

s_t_correct_angle  =  VectorAngle ( s_heading_vector , VectorPointing ( s_geo, t_geo) ) ; 
if  ( s_t_correct_angle  <=  s_fov)  s_encounter  =  true; 

//  Loop  through  rays  and  record  the  ones  seeker  receives 
for  (int  i=0 ; i<d_raycount; i++) { 

//  Initialize  p  variables  to  miss  values 
r [ i ] . power  =  0.0; 

//  Throw  out  rays  that  miss  or  are  obscured 

if  ( r [ i ]  .surf  ==  7  |  |  r[i]  .obscured  ==  true)  continue; 

//**** ***************************GE0MEtry* ************************  ********** 


/*  Check  if  ray  is  within  seeker's  FOV  (does  not  check  reflection  yet) 
The  reason  an  extra  check  is  done  to  see  if  the  ray  is  encounter 
(vs.  target)  is  because  we  may  check  instances  where  the  designator 
greatly  misses  the  target,  and  then  the  seeker  sees  the  overspill/ 
underspill,  even  though  it  is  not  looking  at  the  target  */ 


r_correct_angle  =  VectorAngle ( s_heading_vector, 

VectorPointing (s_geo, r [i ] .hit) ) ; 
if  ( r_correct_angle  >  s_fov)  continue; 

//  Get  Normal  vector  of  surface 
t  normal  vector  =  n[r[i] .surf] ; 

//  Calculate  Distances 

d_t_distance  =  UTILITY . Distance3D (d_geo,  r[i].hit);  // . Desig-target  (m) 

s_t_distance  =  UTILITY . Distance3D ( s_geo,  r[i].hit);  // . Seeker-target  (m) 

//  Calculate  Unit  vectors  Between  Seeker,  Desig,  and  Target  Locations 
s_t_vector  =  VectorPointing ( s_geo,  r[i].hit); 
t_s_vector  =  VectorPointing (r [i] .hit,  s_geo) ; 
d_t_vector  =  VectorPointing (d_geo,  r[i].hit); 

//  Calculate  Angle  between  seeker  and  target  surface  normal 
s_t_normal_angle  =  VectorAngle (t_s_vector,  t_normal_vector) ; 

//  Calculate  Multiplier  for  Ratio  of  Reflected  Radiation  Received 
s_aperture_area  =  M_PI  *  pow ( s_aperture_diameter  /  2,  2); 
s_t_multiplier  =  ( s_aperture_area/ (M_PI  *  pow ( s_t_di stance, 2 )) ) 

*  cos (s_t_normal_angle)  *  cos (r_correct_angle) ; 

//  Break  out  of  loop  if  mult  negative  (surface  normal  >  90  degrees  off) 
if  ( s_t_multiplier  <=  IE-20)  continue; 

//**************************** POWER  LOSS************************************ 
//  Divide  up  total  signal  into  individual  rays 
r[i] .power  =  d_pulse  energy  *  d  efficiency  / 

(d_raycount  *  d_pulse_duration) ; 


if  (r[i]  .surf  ==  6){  // . Ray  goes  directly  from  desig  into  seeker 

double  d_s_distance  =  UTILITY. Distance3D (d_geo,  s_geo) ; 
r[i], power  =  Attenuation (r [i ]. power,  d_s_distance,  d_geo.Z, 

r[i].hit.Z,  sal_weather,  d_lasertype) ; 


} 

else{  // . Ray  reflects  off  of  target  or  ground  before  reaching  seeker 

//  Attenuate  (Designator  to  Target) 

r[i], power  =  Attenuation (r [ i ]. power ,  d_t_distance,  d_geo.Z, 

r[i].hit.Z,  sal_weather,  d_lasertype) ; 

//  Reflect  off  of  Target 

r[i] .power  *=  t_ref lect [ r [ i ] .surf] ; 

//  Attenuate  (Target  to  Seeker) 

r[i], power  =  Attenuation (r [ i ]. power ,  s_t_distance,  r[i].hit.Z, 

s_geo.Z,  sal_weather,  d_lasertype) ; 

//  Decrease  by  distance  and  angle  from  target  normal 
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r[i] .power  *=  s  t  multiplier; 

} 

r[i] .power  *  =  s_ef ficiency; 

//*************DETERMINE  SEEKER  QUADRANT  THAT  RECEIVES  RAY****************** 

//  Define  matrices  to  find  where  ray  intersects  surface 

double  m[3] [3] ; 

double  mi [3] [3] ; 

double  dl[3]; 

double  tuv [3]; 

double  det  =  0.0; 

//  Define  vectors  that  characterize  seeker  pointing 
Point3D  vl  =  s_heading_vector ; 

Point3D  v 2  =  VectorRotation ( s_heading_vector , 0 , 0 , 1 ) ; 

Point3D  v3  =  VectorRotation ( s_heading_vector , 0 , 0 , 2 ) ; 

//  Calculate  m 


m  [  0  ] 

[0] 

=  vl .X; 

m  [  0  ] 

[1] 

=  v2  .  X  ; 

m  [  0  ] 

[2] 

=  v3  .  X 

m  [  1  ] 

[0] 

=  vl.Y; 

m  [  1  ] 

[1] 

=  v2  .  Y  ; 

m  [  1  ] 

[2] 

=  v3.Y 

CM 

£ 

[0] 

=  vl . Z ; 

m  [  2  ] 

[1] 

=  v2 . Z ; 

m  [2  ] 

[2] 

=  v3.Z 

//  Calculate  mi  and  determinant 
InverseMatrix (m, mi , det ) ; 

//  Calculate  dl 


if  (r [i] . surf  ==  6) { 

dl[0]  =  d_geo.X  -  s_geo.X; 
dl [ 1 ]  =  d_geo.Y  -  s_geo.Y; 
dl[2]  =  d_geo.Z  -  s_geo.Z; 

} 


else  { 

dl[0]  =  r[i].hit.X  -  s_geo.X; 

dl[l]  =  r[i].hit.Y  -  s_geo.Y; 

dl[2]  =  r[i].hit.Z  -  s_geo.Z; 

} 

//  Calculate  TUV  vector  (T  =  multiplier,  U  =  pitch,  V  =  yaw) 
tuv [ 0 ]  =  mi [0] [0] *dl [0]  +  mi [ 0 ] [ 1 ] *dl [ 1 ]  +  mi [ 0] [2 ] *dl [2 ] ; 

tuv [ 1 ]  =  mi [1] [0] *dl [0]  +  mi [ 1 ] [ 1 ] *dl [ 1 ]  +  mi [ 1 ] [2 ] *dl [2 ] ; 

tuv [ 2 ]  =  mi [2] [0] *dl [0]  +  mi [2 ] [1] *dl [1]  +  mi [2 ] [2 ] *dl [2 ] ; 

//  Determine  what  quadrant  ray  hits  by  pitch/yaw  corrections 
double  length  =  tuv[0]; 
double  yaw  =  tuv[l]; 
double  pitch  =  tuv[2]; 


//  Sum  up  energies  in  each  quadrant  and  divide  by  pulse  width  to  get  power 
if  (length  >  0) { 
if  (pitch  >  0) { 

if  (yaw  >  0)  {  // . +  Pitch,  +  Yaw 

s^hits [0] ++; 

s_signal.ql  +=  r[i] .power; 

} 

else{  // . +  Pitch,  -  Yaw 

s^hits [1 ] ++; 

s_signal.q2  +=  r[i] .power; 

} 

} 

else  { 

if  (yaw  >  0)  {  // . -  Pitch,  +  Yaw 

s_hits [2 ] ++; 

s_signal.q3  +=  r[i], power; 

} 

else{  // . -  Pitch,  -  Yaw 

s_hits [3] ++; 

s_signal.q4  +=  r[i], power; 

} 

} 

} 
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}/*****************************************************************************/ 
//*******************DETERMINE  DETECT  (S/N  ABOVE  THRESHOLD) ******************* 
double  ql  =  s_signal.ql;  double  q2  =  s_signal.q2; 
double  q3  =  s_signal.q3;  double  q4  =  s_signal.q4; 
double  s_solidangle  =  2  *  M_PI  *  (1  -  cos(sfov)); 

//  Total  seeker  power  is  sum  of  quadrants 
s_power  =  ql  +  q2  +  q3  +  q4; 

//  Calculate  seeker  noise 

s_noise  =  s_background  *  s_aperture_area  *  s_solidangle; 

//  Calculate  each  quadrant's  S/N 
s_signalnoise . ql  =  ql  /  (s^noise  /  4); 
s_signalnoise . q2  =  q2  /  (s^noise  /  4); 
s_signalnoise . q3  =  q3  /  (s^noise  /  4); 
s_signalnoise . q4  =  q4  /  (s_noise  /  4); 

//  Detects  if  at  least  one  quadrant's  S/N  >  threshold 
if  ( s_signalnoise . ql  >  s_signalnoise_threshold  | I 
s_signalnoise . q2  >  s_signalnoise_threshold  | | 
s_signalnoise . q3  >  s_signalnoise_threshold  | | 
s_signalnoise . q4  >  s_signalnoise_threshold) 
s_detect  =  true; 

//***********************  DETERMINE  MANEUVER  S I GNALS  *  ************************** 
if  (spower  <  IE-25) { 

s  maneuver . pitch  =  s  maneuver. yaw  =  0.0; 
smaneuver . theta  =  smaneuver.mag  =  0.0; 

} 

else  { 

s^maneuver .pitch  =  ( (ql  +  q2 )  -  (q3  +  q4 ) )  /  s_power; 
s^maneuver . yaw  =  ( (ql  +  q3)  -  (q2  +  q4))  /  s_power; 

double  p  =  s_maneuver . pitch,  y  =  s  maneuver . yaw; 
int  quadcount  =  0;  double  zerothresh  =  IE-16; 

if  (p  >=  0)  s_maneuver . theta  =  fabs (atan2 (p, -y) ) ; 
else  s_maneuver . theta  =  fabs  (atan2  (p, y) )  +  M_PI; 
if  (ql  <  zerothresh)  quadcount++; 

if  (q2  <  zerothresh)  quadcount++; 

if  (q3  <  zerothresh)  quadcount++; 

if  (q4  <  zerothresh)  quadcount++; 

if  (quadcount  >  1)  smaneuver.mag  =  1.0; 
else  s_maneuver.mag  =  sqrt ( (p*p  +  y*y)  /  2); 

} 

}/*****************************************************************************/ 
/******************************************************************************/ 
/***************★************  UTILITY  FUNCTIONS  *******************************/ 
/********************************************************************* 
double  sSalSeeker Attenuation (double  currentsig,  double  traveldistance, 

double  startheight,  double  endheight, 
weather  sal_weather,  int  lasertype) { 

//  Weather 

int  1  =  sal_weather.  latitude;  // . 0  =  tropics  (<  30  deg  lat) 

// . 1  =  mid-latitudes  (30-60  deg  lat) 

// . 2  =  sub-arctic  (>  60  deg  lat) 

int  s  =  sal  weather . season;  // . 0  =  summer  (3/22-9/21) 
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//.... 

int  v  =  sal_weather . visibility ;  // 
int  la  =  lasertype;  // 

// 

//  Geometry  and  conversion  to  km 
double  td  =  traveldistance  *  .001; 
double  sh  =  startheight  *  .001; 
double  eh  =  endheight  *  .001; 
double  theta  =  asin((eh  -  sh)  /  td)  ; 

//  Attenuation  Coefficient  and  Attenuated  Signal  (Output  of  function) 

double  att  =  0; 

double  attensig  =  currentsig; 


. 1  =  winter  (9/22-3/21) 

0  =  clear  (23  km) ,  1  =  hazy  (5  km) 

. 0  =  Nd:Glass  (1.06um) 

. 1  =  Er:Glass  (1.536um) 


//  Lookup  cables  for  NdtGlass  (1.06  urn)  and  ErtGlass  ( 1 . 536um>  Lasers 
stacic  const  double  L[20] [16]  - 
{ 

//  km~-l 

//  Accen  Coefficient  “  Molecular  Absorption  +  Molecular  Scattering 
//  -t  Aerosol  Absorption  +  Aerosol  Scattering 

//  Rows  0-9  represent  increasing  altitudes  (0  “  0-999m,  9  -  9000-9999m, 
//  heights  >=  10km  set  to  level  9) 


//  Nd:Glass  (1.06 urn.  Note:  Molecular  Absorption  <  le-6,  recorded  as  0) 

/*  - Molecular  Absorption -  - Molecular  Scattering - 

01234S6789  10 


Tr  Summer 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 


Tr  Winter 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 


ML  Summer 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 


ML  Winter 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 


SA  Summer 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 


SA  Winter 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 
0.000000, 


0.000768, 

0.000699, 

0.000633, 

0.000S72, 

0.000519, 

0.000459, 

0.000422, 

0.000380, 

0.000341, 

0.000304, 


0.000768, 

0.000699, 

0.000633, 

0.000572, 

0.000519, 

0.000459, 

0.000422, 

0.000380, 

0.000341, 

0.000304, 


ML  Summer 
0.000781, 
0.000706, 
0.000638, 
0.000577, 
0.000521, 
0.000469, 
0.000421, 
0.000378, 
0.000338, 
0.000302, 


ML  Winter 

0.000843, 

0.000752, 

0.000670, 

0.000599, 

0.000537, 

0.000480, 

0.000427, 

0.000380, 

0.000336, 

0.000297, 


SA  Summer 

0.000798, 

0.000721, 

0.000850, 

0.000584, 

0.000524, 

0.000471, 

0.000423, 

0.000379, 

0.000338, 

0.000301, 


//  Er : Glass  (1.536um) 

/•  - Molecular  Absorption -  - Molecular  Scattering - 

0123456789  10 


Tr  Summer 

0.000131, 

0.000109, 

0.000090, 

0.000075, 

0.000061, 

0.000051, 

0.000041, 

0.000034, 

0.000028, 

0.000022, 


Tr  Winter 

0.000131, 

0.000109, 

0.000090, 

0.000075, 

0.000061, 

0.000051, 

0.000041, 

0.000034, 

0.000028, 

0.000022, 


ML  Summer 
0.000136, 
0.000111, 
0.000092, 
0.000076, 
0.000061, 
0.000051, 
0.000042, 
0.000034, 
0.000026, 
0.000022, 


ML  Winter 

0.000162, 

0.000130, 

0.000104, 

0.000084, 

0.000068, 

0.000055, 

0.000044, 

0.000035, 

0.000028, 

0.000021, 


SA  Summer 
0.000155, 
0.000127, 
0.000104, 
0.000085, 
0.000062, 
0.000057, 
0.000045, 
0.000037, 
0.000030, 
0.000024, 


SA  Winter 

0.000198, 

0.000153, 

0.000120, 

0.000096, 

0.000077, 

0.000062, 

0.000050, 

0.000039, 

0.000030, 

0.000023, 


Tr  Summer 

0.000169, 

0.000154, 

0.000140, 

0.000126, 

0.000114, 

0.000103, 

0.000093, 

0.000084, 

0.000075, 

0.000067, 


Tr  Winter 
0.000169, 
0.000154, 
0.000140, 
0.000126, 
0.000114, 
0.000103, 
0.000093, 
0.000084, 
0.000075, 
0.000067, 


ML  Summer 
0.000172, 
0.000156, 
0.000141, 
0.000127, 
0.000115, 
0.000104, 
0.000093, 
0.000084, 
0.000075, 
0.000067, 


ML  Winter 

0.000186, 

0.000166, 

0.000148, 

0.000132, 

0.000118, 

0.000106, 

0.000094, 

0.000084, 

0.000074, 

0.000066, 


SA  Summer 
0.000176, 
0.000159, 
0.000143, 
0.000129, 
0.000115, 
0.000104, 
0.000093, 
0.000084, 
0.000075, 
0.000066, 


11 

SA  Winter 
0.000877, 
0.000770, 
0.000682, 
0.000606, 
0.000540, 
0.000482, 
0.000429, 
0.000381, 
0.000334, 
0.000288, 


11 

SA  Winter 
0.000193, 
0.000170, 
0.000150, 
0.000134, 
0.000119, 
0.000106, 
0.000095, 
0.000084, 
0.000074, 
0.000063, 


— Aerosol 
12 

Clear 

0.0131, 

0.00571, 

0.00243, 

0.00115, 

0.000723, 

0.000527, 

0.000427, 

0.000418, 

0.000415, 

0.000401, 


— Aercscl 
12 

Clear 

0.0156, 

0.00680, 

0.00290, 

0.00136, 

0.000862, 

0.000623, 

0.000509, 

0.000498, 

0.000495, 

0.000478, 


—Aerosol  Scat - 

14  15 

Clear  Hazy 


Aba— 

13 

Hazy 

0.0582, 

0.0213, 

0.00778, 

0.00284, 

0.00104, 

0.000527, 

0.000427, 

0.000418, 

0.000415, 

0.000401, 


0.0450, 

0.0196, 

0.00836, 

0.00394, 

0.00249, 

0.00181, 

0.00147, 

0.00144, 

0.00143, 

0.00135, 


0.200, 

0.0731, 

0.0267, 

0.00976, 

0.00356, 

0.00181, 

0.00147, 

0.00144, 

0.00143, 

0.00138, 


// 

// 

// 

// 

// 

// 

// 

// 

// 

// 


Abs —  - Aerosol  Scat 

13  14  15 


Hazy 

0.0693, 

0.0253, 

0.00927, 

0.00338, 

0.00124, 

0.000828, 

0.000509, 

0.000498, 

0.000495, 

0.000478, 


Clear  Hazy  »/ 
0.0271,  0.120,  // 
0.0118,  0.0440,  // 
0.00503,  0.0161,  // 
0.00237,  0.00583,  // 
0.00150,  0.00215,  // 
0.00109,  0.00100,  // 
0.000883,  0.000883,// 
0.000865,  0.000865,// 
0.000859,  0.000859,// 
0.000831,  0.000831  // 


0-1  tor. 

1- 2  km 

2- 3  km 

3- 4  km 

4- 5  km 

5- 6  km 

6- 7  km 

7- 8  km 

8- 9  km 

9- 10  km. 


0-1  km 

1- 2  km 

2- 3  km 

3- 4  km 

4- 5  km 

5- 6  km 

6- 7  km 

7- 8  km 

8- 9  tar. 

9- 10  km 


Figure  B-l.  Lookup  tables  for  1.06-  and  1.536-|im  lasers. 

//  Lookup  tables  for  Nd:Glass  (1.06  urn)  and  Er:Glass  (1.536um)  Lasers 
static  const  double  L[20] [16]  = 

{ 

//  kmA-l 

//  Atten  Coefficient  =  Molecular  Absorption  +  Molecular  Scattering 
//  +  Aerosol  Absorption  +  Aerosol  Scattering 

//  Rows  0-9  represent  increasing  altitudes  (0  =  0-999m,  9  =  9000-9999m, 
//  heights  >=  10km  set  to  level  9) 

//  Nd:Glass  (1.06um,  Note:  Molecular  Absorption  <  le-6,  recorded  as  0) 


/*  - Molecular  Absorption -  - 

Molecular  Scattering -  --Aerosol  Abs--  - Aerosol  Scat - 

01234567 
8  9  10  11  12  13  14  15 


Tr  Summer  Tr  Winter  ML  Summer  ML  Winter  SA  Summer  SA  Winter  Tr  Summer  Tr  Winter  ML 
Summer  ML  Winter  SA  Summer  SA  Winter  Clear  Hazy  Clear  Hazy  */ 

0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000768,  0.000768, 
0.000781,  0.000843,  0.000798,  0.000877,  0.0131,  0.0582,  0.0450,  0.200,  //  0-1 

km 

0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000699,  0.000699, 
0.000706,  0.000752,  0.000721,  0.000770,  0.00571,  0.0213,  0.0196,  0.0731,  //  1-2 

km 
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0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000633,  0.000633, 

0.000638,  0.000670,  0.000850,  0.000682,  0.00243,  0.00778,  0.00836,  0.0267,  //  2-3 

km 

0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000572,  0.000572, 

0.000577,  0.000599,  0.000584,  0.000606,  0.00115,  0.00284,  0.00394,  0.00976,  //  3-4 

km 

0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000519,  0.000519, 

0.000521,  0.000537,  0.000524,  0.000540,  0.000723,  0.00104,  0.00249,  0.00356,  //  4-5 

km 

0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000459,  0.000459, 

0.000469,  0.000480,  0.000471,  0.000482,  0.000527,  0.000527,  0.00181,  0.00181,  //  5-6 

km 

0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000422,  0.000422, 

0.000421,  0.000427,  0.000423,  0.000429,  0.000427,  0.000427,  0.00147,  0.00147,  //  6-7 

km 

0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000380,  0.000380, 

0.000378,  0.000380,  0.000379,  0.000381,  0.000418,  0.000418,  0.00144,  0.00144,  //  7-8 

km 

0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000341,  0.000341, 

0.000338,  0.000336,  0.000338,  0.000334,  0.000415,  0.000415,  0.00143,  0.00143,  //  8-9 

km 


0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000000,  0.000304,  0.000304, 
0.000302,  0.000297,  0.000301,  0.000288,  0.000401,  0.000401,  0.00135,  0.00138,  //  9-10 

km 


//  Er:Glass  (1.536um) 


/*  - Molecular  Absorption -  - 

Molecular  Scattering -  --Aerosol  Abs--  - Aerosol  Scat - 

01234567 
8  9  10  11  12  13  14  15 


Tr  Summer  Tr  Winter  ML  Summer  ML  Winter  SA  Summer  SA  Winter  Tr  Summer  Tr  Winter  ML 
Summer  ML  Winter  SA  Summer  SA  Winter  Clear  Hazy  Clear  Hazy  */ 

0.000131,  0.000131,  0.000136,  0.000162,  0.000155,  0.000198,  0.000169,  0.000169, 

0.000172,  0.000186,  0.000176,  0.000193,  0.0156,  0.0693,  0.0271,  0.120,  //  0-1 

km 

0.000109,  0.000109,  0.000111,  0.000130,  0.000127,  0.000153,  0.000154,  0.000154, 

0.000156,  0.000166,  0.000159,  0.000170,  0.00680,  0.0253,  0.0118,  0.0440,  //  1-2 

km 

0.000090,  0.000090,  0.000092,  0.000104,  0.000104,  0.000120,  0.000140,  0.000140, 

0.000141,  0.000148,  0.000143,  0.000150,  0.00290,  0.00927,  0.00503,  0.0161,  //  2-3 

km 

0.000075,  0.000075,  0.000076,  0.000084,  0.000085,  0.000096,  0.000126,  0.000126, 

0.000127,  0.000132,  0.000129,  0.000134,  0.00136,  0.00338,  0.00237,  0.00588,  //  3-4 

km 

0.000061,  0.000061,  0.000061,  0.000068,  0.000062,  0.000077,  0.000114,  0.000114, 

0.000115,  0.000118,  0.000115,  0.000119,  0.000862,  0.00124,  0.00150,  0.00215,  //  4-5 

km 

0.000051,  0.000051,  0.000051,  0.000055,  0.000057,  0.000062,  0.000103,  0.000103, 

0.000104,  0.000106,  0.000104,  0.000106,  0.000623,  0.000828,  0.00109,  0.00100,  //  5-6 

km 

0.000041,  0.000041,  0.000042,  0.000044,  0.000045,  0.000050,  0.000093,  0.000093, 

0.000093,  0.000094,  0.000093,  0.000095,  0.000509,  0.000509,  0.000883,  0.000883,//  6-7 

km 

0.000034,  0.000034,  0.000034,  0.000035,  0.000037,  0.000039,  0.000084,  0.000084, 

0.000084,  0.000084,  0.000084,  0.000084,  0.000498,  0.000498,  0.000865,  0.000865,//  7-8 

km 

0.000028,  0.000028,  0.000026,  0.000028,  0.000030,  0.000030,  0.000075,  0.000075, 

0.000075,  0.000074,  0.000075,  0.000074,  0.000495,  0.000495,  0.000859,  0.000859,//  8-9 

km 
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0.000022,  0.000022,  0.000022,  0.000021,  0.000024,  0.000023,  0.000067,  0.000067, 
0.000067,  0.000066,  0.000066,  0.000063,  0.000478,  0.000478,  0.000831,  0.000831  //  9-10 
km 

}  ; 


//  LEVEL  PATH 
if  (theta  ==  0) { 

//  Figure  out  what  height  laser  is  in 
int  row; 
int  n  =  0; 
do  { 

row  =  n ; 
n+  +  ; 

}  while  (sh  >=  n) ; 
if  (row  >  9)  row  =  9; 


//  Calculate  attenuation  coefficient  for  segment  and  attenuate  the  signal 


//  Molecular  Absorption 

att  =  L[row  +  10*la] [2*1  +  s] 
//  Aerosol  Absorption 

+  L[row  +  10*la] [12  +  v] 
attensig  *  =  exp (-att  *  td)  ; 


Molecular  Scattering 
+  L [ row  +  10*la] [6  +  2*1  +  s] 
Aerosol  Scattering 
+  L[row  +  10*la] [14  +  v] ; 


//  ANGLED  PATH 
else  { 

//  Initialize  new  variables 

double  segment,  ipart;  double  extrasegment  =  0; 
int  n,  lowrow,  highrow; 

//  Make  sure  start  height  is  below  end  height 
if  (theta  <  0) { 

double  sh_temp  =  sh; 
double  eh_temp  =  eh; 
sh  =  eh_temp ; 
eh  =  sh_temp; 

} 


//  Find  altitude  row  that  signal  starts  in 
n  =  0; 
do  { 

lowrow  =  n; 
n+  +  ; 

}  while  (sh  >=  n) ; 
if  (lowrow  >  9)  lowrow  =  9; 


//  Find  altitude  row  that  signal  ends  in 
n  =  0; 
do  { 

highrow  =  n; 
n+  +  ; 

}  while  (eh  >=  n) ; 


if  (highrow  >  9) { 

extrasegment  =  fabs ( (highrow  -  9) /sin (theta) ) ; 
highrow  =  9; 

} 

//  Loop  through  to  attenuate  signal  through  different  altitudes 
for  (int  row  =  lowrow;  row  <  highrow  +  1;  row++) { 

//  Find  out  segment  length 
if  (lowrow  ==  highrow)  segment  =  td; 

else  if  (row  ==  lowrow)  segment  =  fabs((l  -  modf (sh, Sipart) ) /sin  (theta)); 
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else  if  (row  ==  highrow) 

segment  =  fabs (modf (eh, Sipart) /sin (theta) )  t  extrasegment; 
else  segment  =  fabs (1/sin (theta) ) ; 

//  Calculate  attenuation  coefficient  for  segment  and  attenuate  the  signal 
//  Molecular  Absorption  Molecular  Scattering 

att  =  L [ row  +  10*la] [2*1  +  s]  +  L[row  +  10*la] [6  +  2*1  +  s] 

//  Aerosol  Absorption  Aerosol  Scattering 

+  L[row  +  10*la]  [12  tv]  t  L[row  t  10*la]  [14  t  v]  ; 
attensig  *=  exp (-att  *  segment); 

} 

} 

return  attensig; 

j /*****************************************************************************/ 

Point3D  sSalSeeker :: VectorPointing (Point3D  pi,  Point3D  p2,  bool  norm) { 

//  Create  vector  from  pi  to  p2 
Point3D  p3; 

double  vec_mag  =  UTILITY . Distance3D (pi  ,  p2) ; 

//  Subtract  to  get  vector  (p3  =  p2  -  pi) 
p3 . X  =  (p2 . X-pl . X) ; 
p3.Y  =  (p2 . Y-pl . Y) ; 
p3 . Z  =  (p2 . Z-pl . Z) ; 

//  Normalize  vector  if  norm  =  true 
if  (norm  ==  true) { 
p3.X  /=  vec_mag; 
p3.Y  /=  vec_mag; 
p3 . Z  /=  vec_mag; 

} 

return  p3; 

j/*****************************************************************************/ 

Point3D  sSalSeeker : : SpherCartConv (Point2D  pl){ 

//  Creates  unit  vector  from  azimuth/elevation  angles 
Point3D  p3; 

p3.X  =  cos (pi. Y)  *  cos (pi. X); 
p3.Y  =  cos(pl.Y)  *  sin(pl.X); 
p3 . Z  =  sin (pi . Y) ; 
return  p3; 

j/*****************************************************************************/ 

Point3D  sSalSeeker : : VectorRotation (Point3D  pi,  double  a,  double  b,  int  vect) { 
Point3D  p3; 


// . Assume  unit  vectors  for  both  original  and  transformed  vector 

// . pi  =  original  vector,  a  =  yaw,  b  =  pitch 

// . p3_p  =  Perturbed  vector  in  new  prime  coord  system 


double  p3_p[3],  p3_p_mag; 
switch (vect) { 

case  0:  //  Arbitrary  rotations  (default  case) 

p3_p[0]  =  1;  p3_j?[l]  =  tan  (a);  p3_p[2]  =  tan(b);  break; 
case  1:  //  90  degree  rotation  to  ty  axis 

p3_p[0]  =  0;  p3_p[l]  =  1;  p3_p[2]  =  0;  break; 
case  2:  //  90  degree  rotation  to  tz  axis 

p3_p[0]  =  0;  p3_p[l]  =  0;  p3_p[2]  =  1;  break; 

} 

//  Normalize  p3_p [ ] 

p3_p_mag  =  sqrt (p3_p [0] *p3_p [0]  +  p3_p [ 1 ] *p3_p [ 1 ]  t  p3_p [2 ] *p3_p [2 ] ) ; 
p3_p[0]  *=  l/p3_p_mag;  p3_p[l]  *=  l/p3_p_mag;  p3_p[2]  *=  l/p3_p_mag; 

//  Create  rotation  matrix  to  rotate  back  to  original 
//  coordinate  system  using  two  consecutive  Euler  rotations 
//  Rotation  matrix  derived  by  Robert  Yager  in  ARL-TR-5520 
double  R [3 ]  [3]  ; 

double  gamma  =  1/ (sqrt (1-pl . Z*pl . Z) ) ; 
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R [ 0 ] [0]  =  pl.X;  R[0] [1]  =  -gamma *pl . Y;  R[0] [2]  =  -gamma*pl . X*pl . Z ; 

R[l] [0]  =  pl.Y;  R[l] [1]  =  gamma*pl.X;  R [ 1 ] [2]  =  -gamma*pl . Y*pl . Z ; 

R [ 2 ] [ 0 ]  =  pl.Z;  R [ 2 ] [ 1 ]  =  0;  R [2 ] [2]  =  1/gamma; 

p3.X  =  R[0] [0] *p3_p[0]  +  R[0] [1] *p3_p[l]  +  R  [0]  [2]  *p3_p  [2]  ; 
p3 . Y  =  R[l] [0]*p3_p[0]  +  R[l] [l]*p3_p[l]  +  R [1  ]  [2] *p3_p [2]  ; 
p3 . Z  =  R[2]  [ 0 ] *  p  3 _ p [ 0 ]  +  R[2]  [1] *p3_p[l]  +  R [2 ]  [2] *p3_p [2] ; 

return  p3; 

} /**************************************** 
double  sSalSeeker : : VectorAngle (Point3D  pi, 

//  Returns  angle  between  two  vectors 
double  angle; 

//  Calculate  vector  magnitudes  and  dot  product 
double  pl_mag  =  sqrt(pl.X  *  pl.X  +  pl.Y  *  pl.Y  +  pl.Z  *  pl.Z); 
double  p2_mag  =  sqrt(p2.X  *  p2.X  +  p2.Y  *  p2 . Y  +  p2 . Z  *  p2.Z); 
double  dotproduct  =  pl.X  *  p2.X  +  pl.Y  *  p2 . Y  +  pl.Z  *  p2.Z; 
double  division  =  dotproduct/ (pl_mag*p2_mag) ; 

//  Check  boundary  conditions 
if  (division  >  1.0)  division  =  1.0; 
if  (division  <  -1.0)  division  =  -1.0; 


Point3D  p2 ) { 


angle  =  acos (division) ; 
return  angle; 


}/J 


tkkkkkkkkk  < 


double  sSalSeeker :: DotProduct ( Point3D  pi,  Point3D  p2){ 
return  pl.X  *  p2.X  +  pl.Y  *  p2 . Y  +  pl.Z  *  p2 . Z ; 

J.  j  ■k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k  J 

void  sSalSeeker :: Inver seMatrix (double  m[3][3],  double  mi[3][3],  double  &det) { 

/* 

Calculates  inverse  matrix  (mi)  of  matrix  m 
Compute  inverse  matrix  (3x3  specific) 

-1  T 


A  =  label  A  =  1  |  A  B  C  | 

|  d  e  f  |  -  |  D  E  F  | 

I  g  h  i  |  det (A)  | G  H  I  | 


1  I A  D  G  | 

-  |B  E  H| 

det (A) | C  F  I | 


det  =  a(ei-fh)  +  b(fg-id)  +  c(dh-eg) 
*/ 


double  a,  b, 

c. 

d,  e,  f. 

g,  h,  i, 

A,  B, 

c. 

D,  E,  F, 

G,  H,  I; 

a  =  m  [  0  ]  [  0  ]  ; 

b 

=  m  [  0  ]  [1] 

;  c  =  m [ 0 ]  [2]  ; 

d  =  m [  1  ]  [0]  ; 

e 

=  m  [  1  ]  [1] 

;  f  =  m[l] [2] ; 

g  =  m [ 2 ]  [0]  ; 

h 

=  m  [2  ]  [1] 

;  i  =  m [ 2 ]  [2]  ; 

A  =  e*i-f*h; 

B 

=  f*g-d*i 

;  C  =  d*h-e*g; 

D  =  c*h-b*i; 

E 

=  a*i-c*g 

;  F  =  b*g-a*h; 

G  =  b*f-c*e; 

H 

=  c*d-a*f 

;  I  =  a*e-b*d; 

det  =  a*  (e*i 

-f *h)  +  b*  (f 

*g-i*d)  +  c* (d 

mi[0] [0]  =  A/det;  mi[0] [1]  =  D/det;  mi[0] [2]  =  G/det; 

mi [ 1 ] [ 0 ]  =  B/det;  mi[l][l]  =  E/det;  mi[l][2]  =  H/det; 

mi [2 ] [ 0 ]  =  C/det;  mi[2][l]  =  F/det;  mi[2][2]  =  I /det; 

J  j  kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk  J 

void  sSalSeeker :: PrintOutput ( string  filename,  string  format) { 
string  lat,  season,  vis,  laser,  ql ,  q2,  q3,  q4 ,  summary; 
ofstream  outfile; 


//  User  specifies  filename  here  if  not  specified  in  function  call 
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if  (printtofile  ==  true  &&  filename  ==  "") { 
cout  <<  "Please  enter  a  filename:  "; 
cin  »  filename; 
cout  <<  endl  <<  endl; 

filename  =  "OutputW"  +  filename  +  format; 

} 

else  filename  +=  format; 


//  Construct  Quadsignal  Percentages 
if  (s_signal.ql  ==  0.0)  ql  =  ""; 

else  ql  =  ToString (s_signal . ql/s_power  *  100,"%9.2f")  +  "  %" 
if  (s_signal.q2  ==0.0)  q2  =  "  "  ; 

else  q2  =  ToString (s_signal . q2/s_power  *  100,"%9.2f")  +  "  %" 
if  (s_signal.q3  ==  0.0)  q3  =  ""; 

else  q3  =  ToString (s_signal . q3/s_power  *  100,"%9.2f")  +  "  %" 
if  (s_signal.q4  ==  0.0)  q4  =  ""; 

else  q4  =  ToString (s_signal . q4/s_power  *  100,"%9.2f")  +  "  %" 

//  Create  description  strings 
switch (sal_weather . latitude) { 

case  0:  lat  =  "Tropicalin" ;  break; 
case  1:  lat  =  "Mid-lat\n";  break; 
case  2:  lat  =  "Sub-arctic\n" ;  break; 

} 

switch (sal_weather . season) { 

case  0:  season  =  "Summer\n";  break; 
case  1:  season  =  "Winterin'1;  break; 

} 

switch (sal_weather .visibility) { 
case  0:  vis  =  "Clearin'1;  break; 
case  1:  vis  =  "Hazyln";  break; 

} 

switch (d_lasertype) { 

case  0:  laser  =  "Nd: Glassin" ;  break; 
case  1:  laser  =  "Er : Glassin" ;  break; 

} 


//  Construct  Summary  String 

string  fill  =  "  "; 

string  fill2  =  "  "; 

summary  =  "in"  +  Center ("SAL  SEEKER  CLASS  IMPLEMENTATION",  '  ',"") 

+  Center ("developed  for  AWCB  SWEEPM  (LSS;  2011)",  '  ',"")  + 

+  Center ("  INPUT  PARAMETERS  ",  '-',"") 

+  "***  WEATHERin" 

+  "Latitude, 

+  lat 
+  "Season, 

+  season 
+  "Visibility, 

+  vis  +  "in" 

+  "***  SEEKERln" 

+  "Seeker  FOV  (Half  Angle) , 

+  ToString (UTILITY. RAD2DEG  * 

+  "Seeker  Aperture  Diameter, 

+  ToString (s_aperture_diameter  *  1000,  "%-0.2f")  +  "  rnminin 
+  "***  DESIGNATORS" 

+  "Laser  Type, 

+  laser 

+  "Laser  Divergence, 

+  ToString  (d_di-verUence  *  1000, 

+  "Ray  Count, 

+  ToString (d_raycount,  "%-i")  + 

+  "Designator  Error  (H) , 


s_fov,  "%-0.2f")  +  "  degreesin 


"%-0 . 2 f " )  +  "  mradin' 

"in" 
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+  ToString (d_h_error  *  1000,  "%-0.2f")  + 
+  "Designator  Error  (V) , 

+  ToString (d_v_error  *  1000,  "%-0.2f")  + 
+  "***  TARGET\n" 

+  "Size  (Length;  Width;  Height) , 

+  ToString (t_size . X,  "%-0 
+  ToString (t_size . Y,  "%-0 
+  ToString (t_size. Z,  "%-0 
+  "Reflectivity, \n" 

+  "Left  Side, 

+  ToString (t_reflect [0] ,  ' 

+  "Front, 

+  ToString (t_reflect [1 ] ,  ' 

+  "Right  Side, 

+  ToString (t_reflect [2 ] ,  1 
+  "Back, 

+  ToString (t_reflect [3 ] ,  ' 

+  "Top, 

+  ToString (t_reflect  [4 ] ,  1 
+  "Ground, 

+  ToString (t_reflect [5] ,  ' 

+  Center  ("  GEOMETRY  ", 

+  "Designator  to  Target  Distance, 

+  ToString (UTILITY. Distance3D (d_geo, 

+  "Seeker  to  Target  Distance, 

+  ToString (UTILITY. Distance3D( s_geo, 

+  "Seeker  Location  (XYZ  -  m) , 

+  ToString (s_geo .X,  "%-0.0f") 

+  ToString (s_geo.Y,  "%-0.0f") 

+  ToString (s_geo.Z,  "%-0.0f")  +  "\n" 

+  "Designator  Location, 

+  ToString (d_geo .X,  "%-0.0f") 


mrad\n" 

mrad\n\n" 


If' 

')  + 

II 

r 

II 

If' 

')  + 

ii 

r 

II 

If' 

')  + 

"m\n" 

%0. 

.  2  f  " ) 

+ 

”\n 

%0. 

4-) 

CN 

+ 

"\n 

%0. 

.  2  f  " ) 

+ 

"  \  n 

%0. 

.  2  f  "  ) 

+ 

"\n 

%0. 

.  2  f  " ) 

+ 

”\n 

%0. 

.  2  f  "  ) 

+ 

"  \  n 

t_geo) /1000,  "%4.2f")  + 

II 

t_geo) / 1000 ,  "%4.2f")  + 


-O.Of") 


II  O, 

o 

"%-0.0f") 


"%-0.0f") 

"%-0.0f") 

"%-0.0f") 


\n" 


+  ToString (d_geo . Y, 

+  ToString (d_geo . Z, 

+  "Target  Location, 

+  ToString (t_geo .X, 

+  ToString (t_geo . Y, 

+  ToString (t_geo . Z, 

+  "Seeker  Heading, 

+  ToString ( s_heading_vector . X, 

+  ToString ( s_heading_vector . Y, 

+  ToString (s_heading_vector . Z , 

+  "Designator  Heading, 

+  ToString (d_heading_vector . X, 

+  ToString (d_heading_vector . Y, 

+  ToString (d_heading_vector . Z , 

+  "Target  Rotation, 

+  ToString  (UTILITY. RAD2DEG 
+  "Seeker  Azimuth, 

+  ToString (UTILITY. RAD2DEG 
+  "Seeker  Elevation, 

+  ToString (UTILITY. RAD2DEG 
+  Center  ("  SURFACE  HITS  ", 

+  "Surface  Most  Hit, 

+  s__desc  [mainhitsurface]  +  "\n" 

+  "Left  Side, 

+  ToString (surface [0] ,  "%-5i")  +  ", 
+  ToString ( (double) surface [0] *100  / 
+  "Front, 

+  ToString (surface [1 ] ,  "%-5i")  +  ", 
+  ToString ( (double) surface [1 ] *100  / 
+  "Right  Side, 

+  ToString  (surface [2] ,  "%-5i")  +  ", 


\n\n" 

"%-0.2f") 
"%-0.2f") 
"%-0 . 2f " ) 

"%-0.2f") 

"%-0.2f") 

"%-0.2f") 


">\n" 


t_rotation, 
s_orientation . X, 
s_orientation . Y, 


">\n" 
%-0. Of") 


km\n" 

km\n\n" 


<" 


"%-0.0f") 

"%-0.0f") 


degrees\n" 

II 

+  "  degrees\n" 

II 

+  "  degrees\n\n" 


(" 

d_raycount, 

(" 

d_raycount, 

(" 


-O.Of")  +  "% ) \n" 


'%-O.Of")  +  "% ) \n' 
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+  ToString ( (double) surface [2 ] *100  /  d_raycount,  "%-0.0f")  +  "%)\n" 

+  "Back, 

+  ToString (surface [3] ,  "%-5i")  +  ",  (" 

+  ToString ( (double) surface [3] *100  /  d_raycount,  "%-0.0f")  +  "%) \n" 

+  "Top, 

+  ToString (surface [4 ] ,  "%-5i")  +  ",  (" 

+  ToString ( (double) surface [4 ] *100  /  d_raycount,  "%-0.0f")  +  "%)\n" 

+  "Ground,  " 

+  ToString (surface [5] ,  "%-5i")  +  ",  (" 

+  ToString ( (double) surface [5] *100  /  d_raycount,  "%-0.0f")  +  "%)\n" 

+  "DirectlntoSeeker ,  " 

+  ToString (surface [ 6 ] ,  "%-5i")  +  ",  (" 

+  ToString ( (double) surface [6] *100  /  d_raycount,  "%-0.0f")  +  "%)\n" 

+  "Miss, 

+  ToString (surface [7 ] ,  "%-5i")  +  ",  (" 

+  ToString ( (double) surface [7 ] *100  /  d_raycount,  "%-0.0f")  +  "%)\n\n" 
+  Center ("  SEEKER  QUADRANT  HITS  ", 

+  "Q1  (+  Pitch;  +  Yaw) , 

+  ToString ( s_hits [ 0 ] ,  "%-i")  +  "\n" 

+  "Q2  (+  Pitch;  -  Yaw) , 

+  ToString (s_hits [ 1 ] ,  "%-i")  +  "\n" 

+  "Q3  (-  Pitch;  +  Yaw) , 

+  ToString ( s_hits [ 2 ] ,  "%-i")  +  "\n" 

+  "Q4  (-  Pitch;  -  Yaw) , 

+  ToString ( s_hits [ 3 ] ,  "%-i")  +  "\n\n" 

+  Center ("  POWER  LOSS  ",  '-',"") 

+  "Designator  Laser  Signal,  " 

+  ToString (d_pulse_energy  /  d_pulse_duration,  "%-0.2e")  +  "  W\n" 

+  "Seeker  Laser  Signal,  " 

+  ToString (s_power,  "%-0.2e")  +  "  W\n" 

+  "Background  Noise,  " 

+  ToString (s_noise,  "%-0.2e")  +  "  W\n" 

+  "S/N  Threshold,  " 

+  ToString ( s_signalnoise_threshold,  "%-0.2f")  +  "\n" 

+  "S/N  Ql, 

+  ToString ( s_signalnoise . ql ,  "%-0.2f")  +  "\n" 

+  "S/N  Q2 , 

+  ToString ( s_signalnoise . q2 ,  "%-0.2f")  +  "\n" 

+  "S/N  Q3 , 

+  ToString ( s_signalnoise . q3,  "%-0.2f")  +  "\n" 

+  "S/N  Q4 , 

+  ToString ( s_signalnoise . q4 ,  "%-0.2f")  +  "\n" 

+  "S/N  Total7 

+  ToString ( s_power  /  s_noise,  "%-0.2f")  +  "\n\n" 

+  Center ("  ENCOUNTER/ DETECT  ",  '-',"") 

+  "0  =  No;  1  =  Yes\n" 

+  "Target  Encountered,  " 

+  ToString ( s_encounter ,  "%-i")  +  "\n" 

+  "Target  Detected,  " 

+  ToString (s_detect,  "%-i")  +  "\n\n" 

+  Center  ("  QUADRANT  SIGNALS  ",  '-',"")  +  "\n" 

+  fill  +  Center  ("+  Pitch",'  ',"",40) 

+  fill  +  Center  40) ; 

for  (int  n=0;n<3;n++)  summary  +=  fill  +  Center  ("I",'  ’,"|",40); 

summary  +=  fill  +  Align ("Ql  |  Q2",  '  |  ' , 20, " | 40) ; 

summary  +=  fill  +  Center  ("|",'  ',"|",40); 

summary  +=  fill  +  Align (ql  +  "  |  "  +  q2 ,  '  |  ' , 20 , "  |  " , 4 0 )  ; 

for  (int  n=0;n<3;n++)  summary  +=  fill  +  Center  ("|",'  ',"|",40); 

summary  +=  fill2 

+  "+  Yaw  | - |  -  Yaw\n"  ; 

for  (int  n=0;n<3;n++)  summary  +=  fill  +  Center  ("|",'  ',"|",40); 
summary  +=  fill  +  Align ("Q3  I  Q4" , ' | ' , 20, " | " , 40) ; 

summary  +=  fill  +  Center ("|",'  ’,"|",40); 
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summary  +=  fill  +  Align (q3  +  "  |  "  +  q4, ' | , 20, " | " , 40) ; 

for  (int  n=0;n<3;n++)  summary  +=  fill  +  Center  ',"|",40); 

summary  +=  fill  +  Center  4 0) 

+  fill  +  Center  ("-  Pitch",'  ',"",40) 

+  "\n\n" 

+  Center  ("  MANEUVER  GUIDANCE  ",  +  "\n" 

+  "Pitch  Command,  " 

+  ToString ( s_maneuver . pitch,  "%-0.2f")  +  "\n" 

+  "Yaw  Command,  " 

+  ToString ( smaneuver . yaw,  "%-0.2f")  +  "\n" 

+  "Signal  Center:  Theta,  " 

+  ToString (UTILITY. RAD2DEG  *  s_maneuver . theta,  "%-0.2f")  +  "\n" 
+  "Signal  Center:  Magnitude,  " 

+  ToString ( s_maneuver . mag,  "%-0.2f")  +  "\n\n" 

+  Center +  "\n" 


if  (printtoscreen  ==  true)  cout  «  summary  <<  endl; 
if  (printtofile  ==  true) { 

out file . open (filename . c_str ( ) ) ; 
outfile  «  summary  <<  endl; 
out file . close  ( ) ; 

} 

}/*****************************************************************************/ 

main.cpp  #1  (Embeddable) 

/* 

SAL  Seeker  Class  Basic  Implementation  (main.cpp) 

Luke  Strohm,  Army  Research  Laboratory 
3-14-11 
*/ 


/******************************************************************************/ 
#include  <s_sal_seeker_class ,h> 

/******************************************************************************/ 
int  main ( ) 

{ 

sSalSeeker  RUN; 

//****** **************DEFINE  INPUT/OUTPUT  VARIABLES*************************// 

//  INPUT - 

Point3D  d_geo,  t_geo,  s_geo;  // . Designator,  Target,  Seeker 

Point2D  s_orientation;  // . Azimuth,  Elevation  (degrees) 

double  t_rotation;  // . degrees  (0  is  target  traveling  in  the  +y  direct) 

// . rotation  is  clockwise  as  seen  from  above  target 

int  seed;  // . Seed  for  random  number  generator 

//  OUTPUT - 

//  Mega-struct  contains  encounter,  detect,  quad  signals,  and  maneuver  guidance 
saloutputs  results; 

//**** ***************************SET  INPUTS* ************************** ******// 


s 

geo.X  =  -5000; 

s  geo . Y  =  0 ; 

s  geo.Z 

=  1.15; 

d_ 

geo.X  =  -2000; 

d  geo . Y  =  0 ; 

d  geo . Z 

=  1.15; 

t 

geo . X  =  0  ; 

t  geo . Y  =  0 ; 

t  geo.Z 

=  1.15; 

s 

orientation . X 

=  0;  s  orientation . Y  = 

0;  t  rotation  =  0 

seed  =  0; 


//****************************** run  PROGRAM* ********************************// 
//  Convert  inputs  (degrees  to  radians) 
s_orientation . X  *=  RUN . UTILITY . DEG2RAD; 
s  orientation. Y  *=  RUN . UTILITY .DEG2RAD; 
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t_rotation  *=  RUN . UTILITY . DEG2RAD; 

//  Run  main  calculation  function 

results  =  RUN. S_SALSEEKER (d_geo, t_geo, s_geo, s_orientation, t_rotation, seed) ; 

RUN . printtoscreen  =  true; 

RUN . printtof ile  =  false; 

RUN. PrintOutput ("" ,  " . txt") ; 

system ("PAUSE") ; 
return  0; 

} 

main.cpp  #2  (Stand-Alone) 

/* 

SAL  Seeker  Class  implementation  (main.cpp) 

Luke  Strohm,  Army  Research  Laboratory 
3-14-11 
*/ 


/******************************************************************************/ 
#include  <s_sal_seeker_class . h> 

/******************************************************************************/ 
int  main ( ) 

{ 

sSalSeeker  RUN; 

/******************************************************************************/ 
//  Run  Options 
bool  batchrun  =  false; 
bool  testinputs  =  false; 

RUN . printtoscreen  =  true; 

RUN . printtof ile  =  false; 
string  oformat  =  ".txt"; 

/***** *******************DEFINE  input/output  VARIABLES***********************/ 

//INPUT - 

Point3D  d_geo,  t_geo,  s_geo;  // . Designator,  Target,  Seeker 

Point2D  s_orientation;  // . Azimuth,  Elevation  (degrees) 

double  t_rotation;  // . degrees  (0  is  target  traveling  in  the  +y  direct) 

// . rotation  is  clockwise  as  seen  from  above  target 

int  seed;  // . Seed  for  random  number  generator 

//OUTPUT - 

//  Mega-struct  contains  encounter,  detect,  quad  signals,  and  maneuver  guidance 
saloutputs  results; 

/************************* *f>UN  PROGRAM ( S INGLE/BATCH ) *************************/ 
//  Initialize  list  of  file  names  and  find  first  txt  file  in  folder 

_finddata_t  filenamelist;  // . _finddata_t  is  a  Windows-specific  structure 

int  hfile;  // . first  file 

string  ifile,  ofile,  ofile_summary,  outputstring;  // . i/o  files 

of ile_summary  =  "Output\ \OutputSummary . csv" ; 

//  Set  input  folder  and  output  summary  file 

if  (batchrun  ==  true)  hfile  =  _f indf irst ( " InputW* . txt" ,  Sfilenamelist) ; 
else  hfile  =  _findfirst ( "* . txt",  Sfilenamelist) ; 

//  Open  output  summary  file 
ofstream  FileToWrite; 

FileToWrite . open (ofile_summary . c_str ( ) ) ; 
int  counter  =  1; 

do  { 

//SET  UP  I/O  FILES - 
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//  Define  input  and  output  folder  locations 
if  (batchrun  ==  true) { 
ifile  =  "InputW"; 
ifile  =  ifile  +  filenamelist . name; 

} 

else  ifile  =  filenamelist. name; 

//  Open  next  input  file  in  the  folder 
ifstream  FileToRead; 

FileToRead . open (ifile . c_str ( ) ) ; 

//  If  can't  open  input  file,  terminate  with  error 
if  (! FileToRead) { 

cout  <<  "Unable  to  open  file\n"; 
system ("PAUSE") ; 
exit  (1 ) ; 


/**************************QET  DATA  FROM  INPUT  FILE************************/ 
//  Inputs  are  delimited  by  (;) 
for  (int  i=0 ; i<4 4 ; i++) { 

getline (FileToRead,  outputstring,  ';'); 
if  (outputstring  ==  ""){ 

getline (FileToRead,  outputstring) ; 
continue ; 


} 


switch  (i ) 

case 

1  : 

case 

2  : 

case 

3: 

case 

4  : 

case 

5: 

case 

6: 

case 

7  : 

case 

8  : 

case 

9: 

case 

10 

case 

11 

case 

12 

case 

13 

case 

14 

case 

15 

case 

16 

case 

17 

case 

18 

case 

19 

case 

20 

case 

21 

case 

22 

case 

23 

case 

24 

case 

25 

case 

26 

case 

27 

case 

28 

case 

29 

case 

30 

case 

31 

case 

32 

case 

33 

case 

34 

case 

35 

case 

36 

case 

37 

s_geo.X  =  atof (outputstring. c_str ()) ;  break; 
s_geo.Y  =  atof (outputstring. c_str ()) ;  break; 
s_geo.Z  =  atof (outputstring. c_str ()) ;  break; 
d_geo.X  =  atof (outputstring. c_str ()) ;  break; 
d_geo.Y  =  atof (outputstring. c_str ()) ;  break; 
d_geo.Z  =  atof (outputstring. c_str ()) ;  break; 
t_geo.X  =  atof (outputstring. c_str ()) ;  break; 
t_geo.Y  =  atof (outputstring. c_str ()) ;  break; 
t_geo.Z  =  atof (outputstring. c_str ()) ;  break; 
t_rotation  =  atof (outputstring . c_str ()) ;  break; 
s_orientation.X  =  atof (outputstring . c_str ()) ;  break; 
s_orientation . Y  =  atof (outputstring . c_str ()) ;  break; 
seed  =  atoi (outputstring. c_str ()) ;  break; 

break;  // . input  file  line  break 

RUN . sal_weather . latitude  =  atoi (outputstring. c_str ()) ;  break; 
RUN . sal_weather . season  =  atoi (outputstring . c_str ()) ;  break; 

RUN . sal_weather . visibility  =  atoi (outputstring . c_str ()) ;  break; 

break;  // . input  file  line  break 

RUN . s_f ov  =  atof (outputstring . c_str ( ) )  *  RUN . UTILITY . DEG2RAD; 
break; 

RUN . s_aperture_diameter  =  atof (outputstring . c_str ()) ;  break; 

RUN . s_background  =  atof (outputstring . c_str ()) ;  break; 

RUN . s_signalnoise_threshold  =  atof (outputstring . c_str ( ) ) ; break; 
RUN . s_ef f iciency  =  atof (outputstring . c_str ()) ;  break; 

break;  // . input  file  line  break 

RUN . d_lasertype  =  atoi (outputstring . c_str ()) ;  break; 

RUN . d_raycount  =  atoi (outputstring . c_str ()) ;  break; 

RUN . d_divergence  =  atof (outputstring . c_str ()) ;  break; 

RUN . d_pulse_energy  =  atof (outputstring . c_str ()) ;  break; 

RUN . d_pulse_f requency  =  atof (outputstring . c_str ()) ;  break; 

RUN . d_pulse_duration  =  atof (outputstring . c_str ()) ;  break; 

RUN . d_eff iciency  =  atof (outputstring . c_str ()) ;  break; 

RUN . d_h_error  =  atof (outputstring . c_str ()) ;  break; 

RUN . dverror  =  atof (outputstring . c_str ()) ;  break; 

break;  // . input  file  line  break 

RUN.t_size.X  =  atof (outputstring. c_str ()) ;  break; 

RUN.t_size.Y  =  atof (outputstring. c_str ()) ;  break; 

RUN.t_size.Z  =  atof  (outputstring. c  str  ( ) ) ;  break; 
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case  38:  RUN . t_ref lect [0]  = 
case  39:  RUN. t_ref lect [1 ]  = 
case  40:  RUN . t_ref lect [2 ]  = 
case  41:  RUN . t_ref lect [3 ]  = 
case  42:  RUN . t_ref lect [ 4 ]  = 
case  43:  RUN . t_ref lect [ 5 ]  = 

} 

getline (FileToRead 


atof (outputstring. c_str ( ) ) ; 
atof (outputstring. c_str ( ) ) ; 
atof (outputstring. c_str ( ) ) ; 
atof (outputstring. c_str ( ) ) ; 
atof (outputstring. c_str ( ) ) ; 
atof (outputstring. c_str ( ) ) ; 


outputstring) ; 


break; 

break; 

break; 

break; 

break; 

break; 


/**********************DISPLAY  UPDATED  INPUTS  ON  SCREEN********************/ 
if  (testinputs  ==  true) { 


cout 

«  "S  GEO  X 

_  II 

<<  s  geo.X  «  endl 

« 

"S  GEO  Y  =  " 

« 

s  geo.Y  «  endl 

« 

"S  GEO  Z  =  " 

« 

s  geo.Z  «  endl 

« 

"D  GEO  X  =  " 

« 

d  geo.X  «  endl 

« 

"D  GEO  Y  =  " 

« 

d  geo.Y  «  endl 

« 

"D  GEO  Z  =  " 

« 

d  geo.Z  «  endl 

« 

"T  GEO  X  =  " 

« 

t  geo.X  «  endl 

« 

"T  GEO  Y  =  " 

« 

t  geo.Y  «  endl 

« 

"T  GEO  Z  =  " 

« 

t  geo.Z  «  endl 

« 

"T  Rotation 

_  II 

«  t  rotation  «  endl 

« 

"Azimuth  =  " 

« 

s  or ientation . X  «  endl 

« 

"Elevation  = 

II 

<<  s  orientation . Y  «  endl 

« 

"Random  #  Seed 

=  "  <<  seed  <<  endl 

« 

"Latitude  = 

"  «  RUN. sal  weather . latitude  «  endl 

« 

"Season  =  " 

« 

RUN. sal  weather . season  «  endl 

« 

"Visibility 

_  II 

«  RUN. sal  weather . visibility 

«  endl 

« 

"FOV  =  "  « 

RUN 

. S  fov  *  RUN. UTILITY. RAD2DEG  « 

endl 

« 

"Aperture  = 

"  «  RUN.s  aperture  diameter  <<  endl 

« 

"Seeker  Background  Noise  Multiplier  =  "  «  RUN 

.s  background 

« 

"S/N  Threshold 

=  "  <<  RUN.s  signalnoise  threshold  «  endl 

« 

"Seeker  Efficiency  =  "  <<  RUN.s  efficiency  << 

endl 

« 

"Laser  Type 

_  II 

«  RUN.d  lasertype  <<  endl 

« 

"Ray  Count  = 

II 

<<  RUN.d  raycount  «  endl 

« 

"Divergence 

_  II 

«  RUN.d  divergence  «  endl 

« 

"Pulsed  Energy 

=  "  <<  RUN.d  pulse  energy  <<  endl 

« 

"Pulse  Frequency  =  "  <<  RUN.d  pulse  frequency 

«  endl 

« 

"Pulse  Duration 

=  "  «  RUN.d  pulse  duration  « 

endl 

« 

"Designator 

Efficiency  =  "  <<  RUN.d  efficiency 

<<  endl 

« 

"D  H  Error  = 

II 

<<  RUN.d  h  error  «  endl 

« 

"D  V  Error  = 

II 

<<  RUN.d  v  error  «  endl 

« 

"TSize  X  =  " 

« 

RUN.t  size.X  <<  endl 

« 

"TSize  Y  =  " 

« 

RUN.t  size.Y  «  endl 

« 

"TSize  Z  =  " 

« 

RUN.t  size.Z  «  endl 

« 

"TRef lect [ 0 ] 

= 

'  <<  RUN.t  reflect [0]  <<  endl 

« 

"TRef lect [ 1 ] 

= 

'  <<  RUN.t  reflect [1]  <<  endl 

« 

"TRef lect [ 2 ] 

= 

'  «  RUN.t  reflect [2]  «  endl 

« 

"TReflect [3] 

= 

'  «  RUN.t  reflect [3]  «  endl 

« 

"TRef lect [ 4 ] 

= 

'  <<  RUN.t  reflect [4]  <<  endl 

« 

"TReflect [5] 

= 

'  <<  RUN.t  reflect [5]  <<  endl  <<  endl; 

system ("PAUSE") ; 

} 

//  Close  input  file 
FileToRead. close ( ) ; 

//CONVERT  FROM  DEGREES  TO  RADIANS - 

s_orientation . X  *=  RUN . UTILITY . DEG2RAD; 
s_orientation . Y  *=  RUN . UTILITY . DEG2RAD; 
t_rotation  *=  RUN . UTILITY . DEG2RAD; 

/*******************************  CALCULAT IONS*******************************/ 
//  Run  main  calculation  function 
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results  =  RUN. S_SALSEEKER (d_geo, t_geo, s_geo, s_orientation, t_rotation, seed) ; 
/**********************************ouTPUT* ************************ *********/ 


//  Create  output  file  name 

ofile  =  "OutputWSeek  ( "  +  ToString (s_geo .X,  "%-0.0f")  + 

+  ToString ( s_geo . Y,  "%-0.0f")  + 

+  ToString (s_geo.Z,  "%-0.0f")  +  ")" 

+  "Desig("  +  ToString (d_geo . X,  "%-0.0f")  + 

+  ToString (d_geo . Y,  "%-0.0f")  + 

+  ToString (d_geo . Z ,  "%-0.0f")  +  ")" 

+  "Az"  +  ToString ( s_orientation . X  *  RUN. UTILITY. RAD2DEG,  " 

_|_  II  II 

+  "El"  +  ToString ( s_orientation . Y  *  RUN. UTILITY. RAD2DEG,  " 

_|_  II  II 

+  "TR"  +  ToString (t_rotation  *  RUN . UTILITY . RAD2DEG,  "%-0.0 


%-0. Of") 
%-0. Of") 
f")  ; 


//  Output  File  Options 

if  (RUN . pr inttoscreen  ==  true  | |  RUN . printtof ile  == 
RUN. PrintOutput (ofile,  oformat) ; 


true ) 


//  Write  output  to  general  summary  file 
double  totalsig  =  results . sal_signals . ql 
+  results . sal_signals . q2 
+  results . sal_signals . q3 
+  results . sal_signals . q4 ; 
double  totalSN  =  totalsig/RUN. s_noise; 
if  (counter  ==  1 ) { 

FileToWrite  «  "OUTPUT  SUMMARY" 

«  endl  «  endl  «  endl 

«  " - Seeker - ,  ,, - Designator - ,  , 

«  " - Hit  Locations - " 

«  endl 

«  "X,Y,Z,X,Y,Z, Az , El, TR, DHe, DVe, MAIN  HIT," 

«  "FS, F, BS, B, T, G, DirS,Miss, Enc, Det, " 

«  "Q1 , Q2, Q3, Q4, TotSig, TotS/N, Pitch, Yaw, Theta, r" 

«  endl; 

} 

FileToWrite  «  s_geo.X  «  «  s_geo.Y  <<  «  s_geo . Z  «  1 

«  dgeo.X  «  «  dgeo.Y  <<  «  dgeo.Z  «  ' 

«  s_orientation . X  *  RUN . UTILITY . RAD2DEG~« 

«  s_orientation . Y  *  RUN . UTILITY . RAD2DEG  « 

«  t_rotation  *  RUN . UTILITY . RAD2DEG  « 

«  RUN . d_h_error  *  1000  « 

«  RUN . dverror  *  1000  « 

«  RUN . s_desc [RUN . mainhitsur face ]  « 

«  RUN . surface [ 0 ]  «  <<  RUN . surf ace [ 1 ]  << 

«  RUN . surface [2 ]  «  <<  RUN . surf ace [ 3 ]  « 

«  RUN . surface [ 4 ]  «  <<  RUN . surf ace [ 5 ]  « 

«  RUN . surface  [ 6 ]  «  <<  RUN . surf ace [ 7 ]  << 

«  results . encounter  <<  " 

«  results . detect  « 

«  results . sal_signals . ql  *  1000  << 

«  results . sal_signals . q2  *  1000  << 

«  results . sal_signals . q3  *  1000  << 

«  results . sal_signals . q4  *  1000  <<  " 

«  totalsig  *  1000  « 

«  totalSN  « 

«  results . actuator_signals . pitch  « 

«  results . actuator_signals . yaw  « 

«  results . actuator_signals . theta  *  RUN . UTILITY . RAD2DEG  « 
«  results . actuator_signals . mag  « 

«  endl; 


II  II 
r 

ii  ii 
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counter++; 


} 

while  (_f indnext (hf ile, &f ilenamelist)  ==  0)  ;  // . Loop  through  input  files 

//****************************************************************************// 
//  Close  I/O  files 

_findclose  (hfile)  ;  // . Close  the  search  handle 

FileToWrite .  close  ()  ;  // . Close  run  summary  file 

//****************************************************************************// 
if  (RUN . printtoscreen  ==  true)  system ( "PAUSE" ) ; 
return  0; 

}//***************************************************************************// 
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Intentionally  left  blank. 
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Appendix  C.  Utilities  Code 


Utilities  Class 

Written  by  Mary  Arthur  of  the  U.S.  Anny  Research  Laboratory. 

Structs:  Point2D,  Point3D 

Functions:  Distance3D,  DEG2RAD,  RAD2DEG 

Statistics  Class 

Written  by  Mary  Arthur  of  the  U.S.  Anny  Research  Laboratory.  An  implementation  of  the 
routines  found  in 

Saucier,  R.  Computer  Generation  of  Statistical  Distributions,  ARL-TR-2168.  U.S.  Army 
Research  Laboratory:  Aberdeen  Proving  Ground,  MD,  March  2000. 

Functions: 

•  Normal  -  returns  a  random  draw  based  on  a  normal  distribution 

Y  Format  Namespace 

Written  by  Robert  Yager  of  the  U.S.  Army  Research  Laboratory.  Yformat  makes  it  easy  to 
create  an  elegant  output  display,  and  it  was  used  in  the  -PrintOutput’ ’  function  of  sSalSeeker. 

yformatnamespace.h 

J  - k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

: k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 
~k  ~k  ~k  ~k  ~k  ~k 
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**  STRING  FORMATTING  FUNCTIONS  ** 

**  version  1.00  (01-23-2011)  ** 

**  -Rob  Yager  ** 

★  ★  -k  -k  -k  -k 

-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 
•k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k'k-k-k'k-k-k-k-k-k-k-k-k'k-k-k-k'k-k-k'k'k-k-k-k'k-k-k-k-k-k-k  j 

#ifndef  Y_FORMAT_NAMESPACE_H_ 

#def ine  Y_FORMAT_NAMESPACE_H_ 

/******************************************************************************/ 
#include  <string> 
using  std: : string; 

Jki^^ki^-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k/ 

namespace  yFormat{ 

string  Center ( string  s,char  fill='  string  border="  #",int  width=80) ; 
string  Align (string  s, char  center, int  n, string  border="  #",int  width=80) ; 
string  Align2 ( string  s,char  center, int  n, string  s2,int  n2, string  border="  #"); 
string  ToString (double  number, string  format="%f " ) ; 
string  ToString(int  number , string  format="%d" ) ; 
string  ToString ( string  text, string  format="%s" ) ; 

} ;  /  ****************************************************************************/ 
#endif /************************************************************************/ 

y_format_namespace.cc 

#define  _CRT_SECURE_NO_WARNINGS// . disables  deprecation  warnings 

#include  "y_format_namespace . h" 

/ ' k-k'k-k-k'k-k-k-kk'k-k-k'k'k-k-kk'k-k-k'k'k-k-kk'k-k-k'kk-k-k-k'k-k-k'k-k-k-kk'k-k-k'k'k-k-k-k'k-k-k'k-k-k-kk'k-k-k-k-k-k-k-k'k-k-k'k-k-k-kk'k-k-k-k / 

string  yFormat :: Center ( string  s,char  c, string  border, int  width) { 

//reverse  the  order  of  characters  in  border  and  copy  to  rborder 
string  rborder; 

for ( int  i=0; i< (int ) border . size ( ) ; ++i) { 
rborder+=border [border . size ( ) — i — 1 ] ; 

} 

//calculate  the  amount  of  white  space 

double  white_space= ( width- s .size ( ) -2*border . size ( ) ) / 2.0; 

//create  the  centered  string 
string  out; 
out+=border ; 

for (int  i=0 ; i< (int ) white_space ; ++i )  out+=c; 
out+=s ; 

for  (int  i=0;i<(int)  (white_space+0 . 501 ) ; +  +  i )  out+=c; 
out+=rborder+"\n" ; 
return  out; 

j /*****************************************************************************/ 

string  yFormat :: Align ( string  s,char  c, int  n, string  border, int  width) { 
int  count=0; 

for (count=0; count< (int)s.size()&&s [count] ! =c; ++count ) ; //  if ( s [ i] ==c)  count=i ; 
string  out=border; 

for (int  i=0; i<n-count- (int ) border . size () -1 ; ++i)  out+='  '; 
out+=s ; 

int  white spa ce=width- (int) out. size () - (int) border . size ( ) ; 
for (int  i=0 ; i<whitespace ; ++i )  out+='  '; 

//reverse  the  order  of  characters  in  border  and  copy  to  rborder 
string  rborder; 

for ( int  i=0;i<(int)border.size() ; ++i) { 
rborder+=border [border . size ( ) — i — 1 ] ; 

} 

out+=rborder+"\n" ; 
return  out; 

j /*****************************************************************************/ 

string  yFormat : :Align2 (string  s,char  c,int  n, string  s2,int  n2, string  border) { 
int  count=0; 

for (count=0; count< (int)s.size()&&s [count] ! =c; ++count ) ; //  if ( s [ i] ==c)  count=i ; 
string  out=border; 

for (int  i=0 ; i<n-count- ( int ) border . size () -1 ; ++i )  out+='  '; 
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out+=s ; 


int  white spa ce=n2- ( int ) out . size  ( ) ; 
for(int  i=0 ; i<whitespace ; t+i )  out+=' 
outt=s2; 

whitespace=8  0- (int) out .size()-(int) border . size ( )  ; 

//reverse  the  order  of  characters  in  border  and  copy  to  rborder 
string  rborder; 

for ( int  i=0; i< (int ) border . size  ( ) ; ++i) { 
rborder+=border [border . size ( ) — i— 1 ] ; 

} 

outt=rborder+"\n" ; 
return  out; 

}/*****************************************************************************/ 
string  yFormat ToString (double  number , string  format) { 
char  out  [50]  ; 

sprintf (out, format . c_str ( ) , number) ; 
string  string_out=out; 
return  string_out; 

j/*****************************************************************************/ 

string  yFormat ToString ( int  number, string  format) { 
char  out [50] ; 

sprintf (out, format . c_str ( ) , number) ; 
string  string_out=out; 
return  string_out; 

j/*****************************************************************************/ 

string  yFormat ToString ( string  text, string  format) { 
char  out  [200] ; 

sprintf (out, format . c_str ( ) , text . c_str ( ) ) ; 
string  string_out=out; 
return  string_out; 

j/*****************************************************************************/ 
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NO.  OF 

COPIES  ORGANIZATION 


1  DEFENSE  TECHNICAL 
(PDF  INFORMATION  CTR 
only)  DTIC  OCA 

8725  JOHN  J  KINGMAN  RD 
STE  0944 

FORT  BEL  VOIR  VA  22060-6218 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
IMNE  ALC  HRR 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
RDRL  CIO  LL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
RDRL  CIO  MT 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
RDRL  D 

2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 
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NO.  OF 

COPIES  ORGANIZATION 

ABERDEEN  PROVING  GROUND 

37  DIRUSARL 
(35  HC  RDRLWM 
2  CD)  P  PLOSTINS 
M  ZOLTOSKI 
RDRL  WML 
J  NEWILL 
E  SCHMIDT 
T  VONG 
P  WEINACHT 
RDRL  WML  A 
M  ARTHUR 
B  BREECH 
C  MUNSON 
W  OBERLE  (CD  ONLY) 

C  PATTERSON 
R  PEARSON 
L  STROHM  (1  HC,  1  CD) 

A  THOMPSON 
P  WYANT 
R  YAGER 
RDRL  WML  B 
J  MORRIS 
RDRL  WML  C 
BROOS 
RDRL  WML  D 
R  BEYER 
RDRL  WML  E 
I  CELMINS 
F  FRESCONI 
S  SILTON 
RDRL  WML  F 
RHALL 
D  HEPNER 
K  HUBBARD 
P  HUFNAL 
MILG 

G  KATULKA 
D  LYON 
D  PETRICK 
B  TOPPER 
RDRL  WML  G 
T  G  BROWN 
RDRL  WML  H 
B  SORENSEN 
RDRL  WMM 
J  BEATTY 
RDRL  WMP 
P  BAKER 
RDRL  WMS 
T  ROSENBERGER 
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